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Abstract: I give a concise introduction into high energy cosmic ray physics, 
including also few related aspects of high energy gamma-ray and neutrino as- 
trophysics. The main emphasis is placed on astrophysical questions, and the 
level of the presentation is kept basic. 
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1 Introduction 



1.1 Preludes 

What do we want to discuss? 

The term cosmic rays may be defined either as all radiation consisting of relativistic particles 
impinging on the Earth's atmosphere from outer space or, more narrowly, including only 
charged particles. We shall follow generally the latter definition, but discuss also high energy 
gamma rays and neutrinos, since (i) they can be created as secondaries of high energy cosmic 
rays and (ii) they may help us to learn more about the charged component and their sources 
("multi-messenger approach"). The qualification high energy (or relativistic) means that we 
do not consider cosmic rays below a few GeV that may be produced or are infiuenced by the 
Sun and its wind. The cosmic ray spectrum, a nearly featureless power-law extending over 
eleven decades in energy up to a few x 10^'' eV, is shown in Fig. II. 1[ 

We cover mainly astrophysical aspects (What are the sources? How do they accelerate 
cosmic rays? What happens during the journey of cosmic rays to the Earth?), but discuss 
also briefly their interactions in the atmosphere and the experimental methods used to detect 
them. In the final chapter we review connections between cosmic ray physics and searches 
for physics beyond the standard model. 

The level of the presentation is as basic as possible, preferring "back on the envelope 
calculations" , and the material covered corresponds roughly to 20 hours of lectures. 



Where are we? 

Our standard length unit is the parsec, the distance from the Earth to a star that has a 
parallaijll] of one arcsecond. Since one arcsecond is 1/(360 x 60 x 60) = 1/206265 fraction of 
2tt, a parsec corresponds to 206, 265 AU = 3.086 x 10^^ cm = 3.26 lyr = Ipc. 

A schematic picture of our home galaxy, the Milky Way, is shown in Fig. II. 2[ Most stars 
are concentrated in the galactic disc of height h ~ 300 pc in the form of spiral arms. The disc 
is filled with warm atomic gas that consists to 90% of H and to 10% of He and has an average 
density n ~ 1/cm'^. It contains also an ordered magnetic field with strength B ~ 3;uG. The 
energy when the Larmor radius 

R, = ^^ lOOpc ^ ^ (1.1) 

ZeB ^ B Z X IQiSeV ^ ^ 

of a particle with charge Ze and momentum p equals the height of the Galactic disc marks 
approximately the transition between diffusive and rectilinear propagation of cosmic rays. We 
can hope to perform "charged-particle astronomy" only for energies well above this transition. 



^Relatively nearby stars are seen at slightly different positions on the celestial sphere (i.e. the background of 
stars that are "infinitely" far away) as the Earth moves around the Sun. Half of this angular difference is 
called the parallax angle or simply the parallax p. Thus p ~ tanp = AU/d, where 1 AU denotes the mean 
distance between the Earth and the Sun. 
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Figure 1.1: The cosmic ray spectrum I{E) as function of kinetic energy E, compiled using 
results from the LEAP, proton, Akeno, and HiRes experiments jlj. The energy 
region influenced by the Sun is marked yellow and an 1/E'^ '' power-law is also 
shown. 



This is one of the reasons why we are especially interested in cosmic rays of the highest 
energies. 

A halo extends with n ~ 0.01/cm'^ and a turbulent magnetic field probably up to distances 
~ (10 — 15) kpc above the disc. The average strength of this turbulent magnetic field is 
not well restricted and may reach up to < 10 /iG. The whole visible part of the galaxy 
is embedded in a much larger dark matter halo, which comprises 90% of the total mass of 
the Milky Way. The center of the Milky Way contains, as probably most other galaxies, a 
supermassive black hole (SMBH) with mass ~ 10^ Mq. 

The Milky Way is a member of the Local Group, a rather small cluster of galaxies with the 
Andromeda galaxy as the other prominent member. Its diameter is ~ 2Mpc. The distance 
to the nearest galaxy cluster, the Virgo cluster, is 18 Mpc, i.e. approximately equal to the 
mean free path of cosmic rays with the highest observed energies, E ~ lO^'^ eV. 

A naive estimate for the size R of the observable Universe can be obtained from the observed 
recession velocity v = Hod of galaxies at the distance d. Using for the Hubble parameter 
Hq = 70kms~^ Mpc~^, it follows R = c/Hq ~ 4300 Mpc. For comparison, the cosmic ray 
horizon reaches Gpc scale at E = 10^^ eV. 

We shall use mostly natural units with h = c = = ^, but keep often explicitly the speed 
of light c in formulas within a more astrophysical context. 
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Figure 1.2: Schematic picture of the Milky Way with a gas and dust disc, an extended halo 
of gas and cosmic rays, surrounded by globular clusters. Everything is immersed 
in a halo of dark matter. 



1.2 Historical remarks 

1912: Victor Hess discovered on a balloon flight that ionizing radiation increases with altitude. 
As he wrote "The results of the present observations are most easily explained by the 
assumption that radiation with very high penetrating power enters the atmosphere from 
above; even in its lower layers, this radiation produces part of the ionization observed 
in closed vessels. . . Since there was neither a decrease at night or during solar eclipse, 
the Sun can hardly be considered as the source. . . " [2\ . 

1929: Skobelzyn observed first cosmic rays with a cloud chamber. Bothe and Kolhorster 
showed that the tracks are curved by a magnetic field. This proved that the observed 
cosmic rays on ground are charged particles — now we know that these are mainly muons 
produced as secondaries in cosmic ray interactions in the higher atmosphere. 

28/29: Clay observed the "latitude effect": The cosmic ray intensity depends on the (geomag- 
netic) latitude. Bothe and Kolhorster provided first the correct interpretation of this 
effect as an anisotropy induced by the magnetic field of the Earth, providing in turn 
evidence that (the primary) cosmic rays are charged. 

1932: Anderson discovered the positron in cosmic rays. This was just the start for a series 
of new particles detected in cosmic rays: The muon 1936 again by Andersen, charged 
pions in 1947, and strange particles 1947-50. 

1932: Raged debate in the US about sources and primary type of the new radiation. Millikan 
and Compton favored gamma rays and coined therefore the name "cosmic rays." 

1934: The sign of the east- west asymmetry showed that the cosmic ray primaries are positively 
charged particles. 

34/38: Rossi and independently Auger discovered through coincidence measurements "exten- 
sive air showers," showers of secondary particles caused by the collision of high energy 
cosmic rays with air nuclei. 
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1934: Bethe and Heitler developed the electromagnetic cascade theory; the observed particles 

on ground are secondaries. 

1947: Zatsepin discovered the scaling of hadronic interactions studying the evolution of ex- 
tensive air showers. 

1949: Fermi proposed that cosmic rays are accelerated by bouncing off moving magnetic clouds 
in the Galaxy. 

52-54: The first human accelerators reaching p>l GeV were built. As a consequence, cosmic 
ray and high energy physics started to decouple, and cosmic ray physicists focused with 
time more on astrophysical questions. 

1954: First measurements of high energy cosmic rays via sampling of extensive air showers 
done at the Harvard College Observatory. 

1972: The launch of the SAS-2 satellite marked the start of high energy gamma astronomy. 

1976: Start of the first prototype of a large-scale underwater detector for high energy neutrino 
astronomy, DUMAND in Hawaii. 

1998: The Superkamiokande experiment found the first convincing evidence that neutrinos 
are massive observing flavor oscillations of atmospheric neutrinos. 

2007: Completion of the Pierre Auger Observatory in Argentina, the first combination of 
a ground array and fluorescence telescopes in the same experiment. With its size of 
A 3000 km^, it is a factor 30 larger than previous experiments. 
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2 Basic notations of particle physics 



2.1 Relativistic kinematic and cross sections 
2.1.1 Kinematics 

The description of scattering reactions is simplified, if Lorentz invariant quantities arc used. 
Consider for instance the squared center-of-mass energy s = (pa + Pb)^ = {Pc + Pd)^ with 
four- momenta pi = (E'j, pj in the 2^2 scattering a + b ^ c + d, 

s = {Pa+ Pb? = ml + ml + 2{EaEb - PaPb) = ml + ml + 2EaEb{l - PaPb cos i?) , (2.1) 

1/2 

where f3i = Vi = Pi/Ei. We are interested often in the threshold energy s^:^^ of a certain 
process. An important example for cosmic ray physics is the reaction p + 7 — p + tt^. If in 
this process the proton is at rest, then 

s = + 2E^mp > {rup + m.,^)'^ = + 2mpmT^ + (2.2) 

or 

2 

E^>m^ + « 145 MeV . (2.3) 

2mp 

Thus the photo-production of pions on protons at rest is only possible for E^ > 145 MeV. 
Photons with lower energy do not interact via this reaction, while higher energetic ones do. 
We will consider this process later again, but then as the scattering of ultrahigh energy protons 
on low energy photons, that is an important example for the transition from a transparent 
Universe at low energies to an opaque one at high energies. 

Ex. : Calculate the minimal energy of a proton able to produce anti-protons scattering on another 
proton at rest. 

Because of baryon number conservation, anti-proton production is first possible in pp pppp with 
s = 2nip + 2Epmp > IGnip or Ep > Tnip. Furthermore, the cross section of this reaction is small 
close to the threshold. Hence the anti-proton flux at small energies should be strongly suppressed, if 
anti-protons arc only produced as secondaries in cosmic ray interactions. 

The second important quantity characterizing a scattering process is the four-momentum 
transfer t = {pa — Pc)^ = (Pb — Pdf- As an example, we consider electron-proton scattering 
e~ + p ^ e~ + p. Then 

t={Pe- P'e? = - 2Ee4 (1 - (iefi'e COS 1») . (2.4) 

For high energies, Pe^Pe ~^ 1; ^i^d 

t ?s -2EeE'^{l - cos -d) = -4.EeE'^ sin^ 'd/2 . (2.5) 

Because of four-momentum conservation, the variable t corresponds to the squared momentum 
Q of the exchanged virtual photon, t = = {pe — p'^,)^ < 0. Thus a virtual particle does 
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not fulfill the relativistic energy-momentum relation. The energy-time uncertainty relation 
AEAt > 1 allows such a violation, if the virtual particle is exchanged only during a short 
enough time. Note also that the angular dependence of Rutherford scattering, dcj/dO oc 
l/sin^i?/2 is obtained if da/dQ cx l/t^. 

In the reaction e~ + p ^ e~ + p the scattering happens most likely at small angles ■(? and 
small four-momentum transfer t. (Note that this does not imply small energy transfer in the 
lab system.) In a collider experiment, this kinematical region is difficult to observe because 
of the beam pipe. Moreover, the search for physics beyond the standard model as one of 
the most important problems in high energy physics is often performed as search for new, 
heavy particles. The production of particles with mass M requires typically \t\ > and 
is therefore rare. On the other hand, for the description of cosmic ray interactions in the 
atmosphere we are interested mainly in the most common interactions. This is one of the 
reasons why high energy particle physics with accelerators and with cosmic rays are focusing 
partly on different aspects of high energy interactions. 



2.1.2 Cross section, interaction depth, and slant depth 

Consider a cylinder of length I and area A filled with N scattering centers each with cross 
section a. Their number density n is n = N/{Al). Let us assume first that Na <^ A. Then 
the fraction of incoming particles absorbed in the cylinder is simply 

Na , , 

— = nla = T (2.6) 

and defines the optical or interaction depth r = nla. Our assumption Na <C A corresponds 
to T <C 1, an "(optical) thin" or transparent source in contrast to an "(optical) thick" source 
with T » 1. In general, we have to calculate how much radiation is absorbed passing the 
infinitesimal distance d/, 

dl = -IdT = -Inadl . (2.7) 

Integrating gives 

I = Io exp(-r) or = Iq exp ^- ^ dl na^ . (2.8) 

In cosmic ray physics, one introduces often the interaction length Aj = m/ai measured in 
g/cm'^ and replaces the path length I by the slant depth X, 

X{1)= / dlp{l') or Xih)= / dhp{h'). (2.9) 

JO Jh 

Hence, X measures the weight per area of the material crossed. The version on the RHS of 
Eq. (|2.9p is appropriate in the case of the Earth atmosphere, where X = corresponds to 
h — > oo, while X measured from the zenith is numerically equal to 1030 g/cm^ at sea level. 
The advantage of X is obviously to hide the integration in X, 

A A 
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Table 2.1: The known particles around 1950 together with their main decay mode, their 
life-time and range. 

2.2 Particles and interactions 

2.2.1 Our particle inventory 

The particles known around 1950 are shown in Table [2?T] together with their main decay mode 
and their life-time ti/2 and range cti/2- Some immediate consequences of the life-times given 
are: 

• Neutrons with T = E/nin > 10^ or E > 10^^ eV are stable on galactic scales (~ lOkpc). 
A source of neutrons with E > 10^^ eV in our galaxy would provide neutral, strong 
interacting primaries, while at lower energies it would be visible as a i>e source from 
neutron decay on the flight. 

• In the atmosphere, only low energy charged pions will decay and produce neutrinos, 
whilst at high energies charged pions mainly scatter. Similarly, most high energy muons 
do not decay within the extension of the atmosphere, ~ 15 km. Thus the neutrino flux 
produced by cosmic rays in the atmosphere should be a steeper function of energy 
than the cosmic ray flux. This opens the possibility to perform high energy neutrino 
astronomy at energies E ^ TeV, for which the background of atmospheric neutrinos 
becomes negligible. 

• The long muon range can be used to enlarge the effective volume of a neutrino detector, 
observing also muons produced outside the proper detector volume. 

2.2.2 Comparison of electromagnetic, weak and strong interactions 

The energy-time uncertainty relation AEAt > 1 restricts the emission of a heavy particle 
with mass M by a lighter one to times At < l/M. Thus the range At of an interaction over 
which a massive particle can be exchanged should be limited by l/M. 

This idea leads to the generalization of the Coulomb potential to the Yukawa potential, 

V{r) = ^"P(-^^) ^ I for M^O. (2.11) 

In non-relativistic quantum mechanics, the connection between the cross section a and the 
potential V is given via the scattering amplitude f{'&), 

^ = l^^^ = l d^^l/Wl'- (2-12) 
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In the so-called Born approximation, this scattering amplitude /(i?) is the Fourier transform 
of the potential, 

f{Q) = / d'l exp(-z(qr))F(r) , (2.13) 

where q = p — p' or = 2psin'!?/2. 

Performing the integral for the Yukawa potential one obtains 

dn (g2 + M2)2 ■ ^^-^^^ 

This result explains the difference in strength between electromagnetic and weak interactions, 
although their coupling constants have roughly the same size, 5 ~ e. Photons are massless 
and thus 

while the bosons exchanged in weak interactions are heavy. If s ^ M^, then also |t| ^ 
and ^ M^, and thus the exchange of massive gauge bosons results in 

^«^«G|, w.th G..^. (2.16) 

The masses of the W and Z bosons that intermediate weak interactions have been measured 
at accelerators as Mz = 91.2 GeV and Mw = 80.4 GeV. 

The total cross section a is obtained by integrating over d$7 oc d(cosT?) oc dt. We expect 
from dimensional analysis that a ~ G|n (energy)^ and often the result can be estimated. For 
instance, for the scattering of a high-energy neutrino on a nucleon at rest the integration gives 
an additional factor s 2mjv-Ei/, and the cross section is numerically for s ^ 

.„.~OJ»»10-3W(^^). (2.17) 

Finally, we consider strong interactions. The largest contribution to the proton-proton 
cross section should result from the exchange of the lightest, strongly interacting particle, the 
pion. Integrating gives 



g4 /-O g Q,2 

a ~ / 7: ^ = al ^ ^ ^ « 200^ ^^1^^^,^^ (2.I8) 

J-s (t — m^y m^{s — m^j 



(47r)2 {t - m2)2 ^ m|(s - ml) 

compared to the measured value a ^ 50mbarn. Hence the proton-pion coupling = 5^/(47r) 
is of order one and perturbation theory for strong interactions seems to make no sense. 
However, coupling constants are — despite the name — in generally not constant but depend 
on the scale probed. In the case of strong interactions, the coupling becomes smaller 
for large Q^, a phenomenon called "asymptotic freedom." As a consequence, perturbation 
theory can be used in QCD, however only in a restricted kinematical range, > few x GeV^. 
Unfortunately, for cosmic ray interactions most important is the regime of small momentum- 
transfer that is not (yet) accessible to predictions from first principles and requires therefore 
the use of phenomenological models. 

We have argued that the exchange of the lightest (relevant) particle dominates the total 
cross section. But what happens, when an infinite number of particles can be transmitted. 
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Table 2.2: Valence quark content of hadrons 



with increasing multiplicity for increasing mass? This is the case of hadrons for which a zoo 
of resonances exists, cf. the Chew-Prautschi plot in Fig. 12.11 For instance, the vector mesons 
A2, ps, ■ ■ ■ can be viewed as J = 2, 3, . . . excitations of the p meson with spin J = 1. After 
summing up the infinite number of exchanged hadrons per trajectory, one can replace them 
by a single "effective particle" (called "reggeons"). The high energy behavior of hadronic 
cross sections is modeled as an exchange of these reggeons, with their properties read off from 
Chew-Frautschi plots as in Fig. 12. 1[ 



Quark model and hadronization For q > l/nip, the scattering amplitude for e.g. electron- 
proton scattering should reflect the proton size 10~^^ cm. However, a form factor typical for 
scattering on point-like particles was found experimentally. This led to the parton model: 
Protons (as all hadrons) consist of point-like constituents, quarks and gluons. The valence 
quark content of some hadrons is given in Tab. 12. 2i At high momentum transfer, a probe 
scatters on partons independently ("spectator model"). Hence, the parton model predicts 
e.g. for the ratio of pion-proton and proton-proton cross sections 

^ = I , (2.19) 

in good agreement with experiment. 

Which kind of hadrons are mainly produced in the process e~^e~ — > qq —> hadrons? Mainly 
mesons (qq) with only ~ 5% of nucleons (qqq), because it is easier to combine two quarks into 
a colorless, light meson than three quarks into a colorless, rather heavy baryon. Out of the 
mesons, 90% are pious while the remaining ones are mainly kaons. While 1/3 of pions are 
neutral ones and produce photons, the other 2/3 are charged and produce neutrinos. Hence 
photon and neutrino fluxes produced by hadronization are closely connected. 



2.3 Exercises 

1. Derive the connection between da/dO, and da/dt and the integration limits in t for a 
general 2 — > 2 scattering process. 

2. Find the minimal energy E'th of a proton scattering on a photon with the typical energy 
of the cosmic microwave background (T « 2.7 K) for the process p + j^p + ir^. Guess 



the cross section of this reaction, check it against the curve at http://pdg.lbl.gov 
and estimate the mean free path of a proton with E ^ E'th. 

Consider the 2-particle decay vr*^ — > 27. What are the minimal and maximal photon 
energies, if the pion moves with velocity v? What is the shape of the photon spectrum 
diV/dE? 
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Figure 2.1: Chew-Frautschi plot, i.e. the spin against the squared mass, for the trajectories 
of the w, p, f2, and a2 mesons, from Ref. [3j. 
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3 Galactic cosmic rays 



3.1 Basic observations 



3.1.1 Cosmic ray intensity and composition 

The integral intensity /(> E) of cosmic rays is defined as the number of particles with 
energy > E crossing an unit area per unit time and unit solid angle. Thus its units are 
[/] = cm~^s~^sr^^. The (differential) intensity I{E) and the integral intensity are connected 
by 

roo 

/(> E) = / dE'l{E'). (3.1) 

JE 

The particle flux from one hemisphere through a planar detector is 



27r 



T{E) = / dn I{E) cos 1? = I{E) d(t> d'& sin i? cos i? = 7rI{E) 



■k/2 



tt/2 



di? sin 2i? 



7Tl{E) , 

(3.2) 

where we have assumed that the intensity is isotropic. The (differential) number density of 
cosmic rays with velocity v is 

n{E) 



47r 

-m . 

V 



(3.3) 



More generally, the intensity I may depend both 
on the position x of the detector and on its ori- 
entation (p, We can connect / to the phase 
space distribution /(x, p) of cosmic rays by com- 
paring dA^ = /(x, p)d^xd^p with the definition 
/ = dN/{dA dt dl^dE), 



/(x,p,??,, 



2 dp . 
vp — /(x,p) 



pV(x,p). (3.4) 



Figure [3T] shows I{E) separately for some im- 
portant elements. First, one recognizes that the 
main component of cosmic rays are protons, with 
additionally around 10% of helium and an even 
smaller admixture of heavier elements. 

The relative abundance of elements measured 
in cosmic rays (dark, filled circles) is compared 
to the one in solar system (blue, open circles) in 
Fig. 13.21 Both curves show the odd even effect, 
i.e. the tighter bound nuclei with an even numbers 
of protons and neutrons are more abundant. The 
main difference of the two curves is that the Li- 
Be-B group (Z = 3 - 5) and the Sc-Ti-V-Cr-Mn 
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Figure 3.1: Cosmic rays intensity for sev- 
eral elements, from Ref. [1]. 
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[Z = 21 — 25) group are much more abundant in 
cosmic rays than in the solar system. We explain 
this later as a propagation effect: The elements 

from the Li-Be-B group are produced as secondaries in cosmic rays interactions in the Galaxy, 

Second, the spectra shown in Fig. 13. H are above 
a few GeV power-laws, practically without any 
spectral features. The total cosmic ray spectrum 
is 

particles 



I{E) ~ l.SE' 



cm^ s st GeV 



(3.5) 



in the energy range from a few GeV to 100 TeV 
with a « 2.7. Around 10^^ eV (the "knee"), the 
slope steepens from a « 2.7 to a « 3.0. The spec- 
trum above 10^^ eV will be discussed in Chap- 
ter 6. 

The power-law form of the cosmic ray spec- 
trum indicates that they are produced via non- 
thermal processes, in contrast to all other radi- 
ation sources like e.g. stars or (super-) novae 
known until the 1950's. 

Third, small differences in the exponent a of 
the power-law for different elements are visible: 
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Figure 3.2: Abundance of elements mea- 
sured in cosmic rays compared 
to the solar system abun- 



dance, from Ref. [5]. 
The relative contribution of heavy elements increases with energy. 
The kinetic energy density of cosmic rays is 



PCR = JdE Ekn{E) =47TjdE-^ I{E) 



(3.6) 



Extrapolated outside the reach of the solar wind, it is 



PGR ~ 0.8eV/cm^ 



(3.7) 



compared to the average energy density pb ^ lOOeV/cm'^ of baryons in the Universe, of star 
light /Oiight ~ 5eV/cm^ in the disc, and in magnetic fields pmag = 0.5eV/cm^ for B = 6/iG. 
If the local value of pcR would be representative for the Universe, 1% of the energy of all 
baryons would be in the form of relativistic particles. This is rather unlikely and suggests 
that cosmic rays are accumulated in the Galaxy. 



Solar modulations When cosmic rays enter our Solar System, they must overcome the 
outward- flowing solar wind. This wind impedes and slows the incoming cosmic rays, reducing 
their energy and preventing the lowest energy ones from reaching the Earth. This effect is 
known as solar modulation. The Sun has an 11-year activity cycle which is reflected in the 
ability of the solar wind to modulate cosmic rays. As a result, the cosmic ray intensity at 
Earth is anti-correlated with the level of solar activity, i.e., when solar activity is high and 
there are lots of sunspots, the cosmic ray intensity at Earth is low, and vice versa. 

Since the number of cosmic rays increases with decreasing energy, most cosmic rays are not 
visible to us. This suppression effect at energies below a few GeV is clearly visible in Fig. 13.31 
where the intensity of oxygen is shown for three different periods. 
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Figure 3.3: Oxygen cosmic ray intensity during three different periods: Sept. 1997 (squares), 
Feb 2000 (circles), Jan. 2001 (diamonds), from Ref. [6]. 



3.1.2 Anisotropies and deflections in regular magnetic fields 

Experimental results on dipole anisotropies 

Above 100 GeV, when the influence of the Sun 
becomes negligible, the cosmic ray flux is consis- 
tent with isotropy. At higher energies, when cos- 
mic rays drift slowly out the galaxy, one expects 
to detect first anisotropies on large angular scales. 
Such anisotropies can be characterized in first ap- 
proximation by a pure dipole, I = Iq + Ii cost?, 



-§10 

"q. 
E 
< 



o 


Ml. Norikum 1990 


□ 


EAS -TOP 1996 


A 


EAS-TOP 2003 




Akeno 1986 




Adelaide 19S3 



10 ^ 



KASCADE 



In 



+ In 



h 



(3.8) 

Most experimental detections or limits for 5 
at energies E < 10^^ eV are in the range 
10~^ < S < 10~^, cf. Fig. 13. 4[ Many experimen- 
tal searches are using that the experimental ex- 
posure of a cosmic ray experiment taking contin- 
uously data is uniform in right ascensioru. Then 
one performs an one-dimensional harmonic anal- 
ysis: One sums 

2 " 2 
Ofc = — cos(k(j)i) and = — sin(k(f)i) 
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Figure 3.4: Experimental limits on the 
anisotropy of cosmic rays, 
from Ref. 13 • 



n 



n 



a=l 



(3.9) 



^In the equatorial coordinate system one projects the usual geographical latitude and longitude on the celestial 
sphere. The right ascension measures the angle east from the vernal equinox point to an object projected 
on the celestial equator. 
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over the n data points. The amphtude and and phase (pk of the k.th harmonic are given 

by 



rk = Y^fc + bl and (pk = arctan(6fc/afc) (3.10) 
with the chance probabiUty 

Pch = exp (^-nrl/4j (3.11) 
to observe a larger value of in an isotropic distribution. 

Generalized Liouville theorem We measure cosmic rays after they traversed the magnetic 
field of the Earth and of the Milky Way, and possibly extragalactic magnetic fields. It is 
therefore important to separate between genuine anisotropics and those eventually induced 
by the cosmic ray propagation in magnetic fields. 

We want to show now that the intensity / is constant along any possible cosmic ray tra- 
jectory. Let us consider the evolution of the phase space distribution /(x, p) from the time t 
to t + dt. The number dA^ of particles around the point x', p' at t' = t + dt is 

/(x',p')d3x'dV. (3.12) 
This number remains constant, if the Jacobian of the transformation x, p ^ x',p' is one, 

= (3.13) 

9(x,p) 

It is sufficient to show that the time-derivative dJ/dt = 0, i.e. that in the expansion of J the 
first-order terms dt vanish. From x' = x + vdt and p' = p + Fdt, the diagonal of J follows as 

8F BF BF 

1, 1, 1, 1 + ^dt, 1 + pdt, 1 + ^di, (3.14) 

OPx Bpy dpz 

while the off-diagonal elements are of order 0{dt). The expansion of dJ to first-order in dt 
is thus 

f BF^ BFy BFA , , 

y Bpx Bpy Bpz J 

Thus the phase space distribution /(x, p) is constant along a trajectory, if VpF = 0, which 
is fulfilled in an electromagnetic field. Moreover, p = |p| is constant in a pure magnetic field 
and hence also p^ f{^,p). With Eq. (j3.4p it follows that also the intensity / is constant along 
any possible cosmic ray trajectory. 

Finally, we remark that the Liouville theorem does not imply that magnetic fields can 
not enhance existing anisotropics. In fact, they can act as magnetic lenses, enhancing or 
decreasing the fiux received from point sources similar as gravitational lenses. 



Latitude and east-west effect An isotropic cosmic ray fiux remains isotropic propagating 
through a magnetic field as long as the phase space is simply connected. In other words, a 
necessary condition is that all trajectories starting from the point considered on Earth (after 
reversing the charge of the particle) reach r = oo. At low enough energies, this condition 
may be violated, because trajectories can be defiected back to the Earth or stay within a 
finite distance r, cf. Fig. 13. 5i In this case, the magnetic field does induce anisotropics in the 
observed fiux. 
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Figure 3.5: A sketch of the east-west effect [S]. 



Consider a particle of charge Ze with orbit in the equatorial plane of a dipole with magnetic 
moment M (as a good model for the geomagnetic field). Equating the centrifugal and the 
Lorentz force gives 

2 

TTIV 

Ze\v X B| = (3.16) 



r 



with B = /io/(47r)M/r''. The radius of the orbit is 



f«l^V' . (3.17) 

V47r p J 



Setting r = -R® and using M = 8 x 10 Am as magnetic moment of the Earth, it follows 

This is the minimal momentum of a proton able to reach the Earth from the east, if its orbit 
is exactly in the (magnetic) equatorial plane. The sign of this east-west asymmetry was used 
by Rossi and others to show that the cosmic ray primaries are positively charged. Towards 
the poles, the influence of the dipole field becomes weaker (v x B), and the cutoff momentum 
becomes thus smaller. Thus the integrated cosmic ray intensity increases with latitude for 
charged particles ("latitude effect"). 



Compton-Getting effect Compton and Getting first discussed that a relative motion of 
observer and cosmic ray sources results in an anisotropic cosmic ray flux, using this effect as 
signature for the Galactic origin of cosmic rays with > 0.1 GeV. 

Lorentz invarianc^ requires that the phase space distribution function / in the frame of 
the observer, /'(r', p'), equals the one in the frame in which the cosmic ray flux is isotropic, 

■^The differentials d^x and d^p transform opposite under Lorentz transformations, while the particle number 
dA'^ is obviously a scalar. 
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/(r, p). Expanding in the small parameter p — p' ~ — pu, it follows 

,V)^/,p')-.u.MEO..(„^,./,p')(l-^ii^)^ (3.., 

Since u= |u| <C 1, the anisotropy induced by the Compton-Getting effect is dominated by the 
lowest moment, i.e. its dipole moment. Changing to the differential intensity I{E) = p'^f{p), 
one obtains 

d In / \ u • p " 



,'(E')^m\l+[2-^):^y (3.20) 

Thus the dipole anisotropy due to the Compton-Getting effect has the amplitude 

OCG = J —j: — = 2-— ^ -u. (3.21) 

-'max + -imin V Q m / 

The Sun moves with uq = 220km/s around the center of the Milky Way. Most likely, the 
local "cosmic ray rest frame" is co-rotating with the nearby stars and the relevant velocity 
u for the CG effect is therefore much smaller. Taking into account the observed spectrum 
1(E) oc E~'^''^ of cosmic rays below the knee, the Compton-Getting effect should results in a 
dipole anisotropy which amplitude is bounded by 5cg = (2 -|- 2.7) u < 0.4%. 

3.2 Propagation of Galactic cosmic rays 

Cascade equation We want to explain the large over-abundance of the group Li-Be-B in 
cosmic rays compared to the Solar system. We consider two species, primaries with number 
density Up and secondaries with number density n^. If the two species are coupled by the 
spallation process p ^ s + X, then 



dup Up 



(3.22a) 



where X = J dl p{l) measures the amount of traversed matter, Aj = m/ai are the interaction 
lengths (in gr/cm^), and p^p = Cgp/cxtot is the spallation probability. 

The system Eqs. (j3.22p can be easily solved (Exercise 3.1) and using as initial condition 
n<j(0) = we obtain as ratio 



As — Ap 




exp T 1- - 1 



(3.23) 



If we consider as secondaries a group like Li-Be-B that has a much smaller abundance in 
the solar system than in cosmic rays, most of them have to be produced by spallation from 
heavier elements like the C-N-0 group. With Acno ~ 6.7g/cm^, ALiBeB ~ lOg/cm^, and 
Psp ~ 0.35 measured at accelerators, the observed ratio 0.25 is reproduced for X ~ 4.3g/cm^, 
see Fig. 13.61 

With h = 300 pc ~ 10^^ cm as thickness of the Galactic disc, uh ~ 1/cm'^ as density of 
the interstellar medium, a cosmic ray following a straight line perpendicular the disc crosses 
only X = mnn-Hh ~ 10~'^g/cm^. The residence time of cosmic rays in the galaxy follows 
as t ~ (4.3/10~^)(/i/c) ~ 1.4 x lO^^s ~ 5 x 10^ yr. This result can only be explained, if the 
propagation of cosmic rays resembles a random- walk. Moreover, it suggests that acceleration 
and propagation can be treated separately. 



20 



3.2 Propagation of Galactic cosmic rays 




Figure 3.6: The ratio rig /rip as function of the traversed amount of matter X. 

Random walks After N steps li of the same size |li| = / a particle that started at zero is at 
the position d = J2iLi !«• We assume that the direction of each step li is chosen randomly. 
Then the scalar product of d with itself is 

N N 

d.d = J2J2^,.l,, (3.24) 

i=ij=i 

and splitting the sum into the diagonal and the off-diagonal terms, we obtain 

N N 

dP = + 2f^^cosdij ^ Nl'^ . (3.25) 

1=1 j<i 

By assumption, the angles "dij between Ij and Ij are chosen randomly and thus the off-diagonal 
terms cancel against each other. 

One might think of two different reasons behind the random walk: Scattering with isotropic 
re-emission as in the case of photons in the Sun or diffusion in turbulent magnetic fields. The 
first possibility is excluded by the small scattering probability r = hann ~ 10^^ x 10~^^ x 1 ~ 
10~^ and by the fact that at high energies scattered particle are concentrated in the forward 
direction within a cone of opening angle ~ l/F. Hence we are lead to examine the second 
possibility, the scattering of cosmic rays on turbulent magnetic fields in the Galactic disc. 



DifFusion equation The continuity equation for the number density n and its current j (that 
corresponds to JT in our usual notation), 

„. dn 

Vj + — = 0, (3.26) 
leads together with Tick's law for an isotropic medium, 

j = -DVn, (3.27) 

to the diffusion equation 

UTl 

— -V{DVn) = Q, (3.28) 
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where we added additionally a source term Q = Q{E,x,t). Note that the diffusion equa- 
tion can be transformed for a diffusion coefficient D that is independent on x into the free 
Schrodinger equation substituting D <-> fi^ /{2m) and t <-> —it. Hence we can borrow the free 
propagator for a non-relativistic particle as Green's function G{r) for the diffusion equation 
with D = const, and obtain with the mentioned substitutions, 

Gir) = j^:;^,^M-rV{m\- (3-29) 

Thus the mean distance traveled outward is oc V Dt, as in a random walk with (r^) ~ -^^o- 
Connecting the two pictures, we obtain D ~ Nl^/t ~ vIq with v = Nl^/t. Therefore, the 
diffusion coefficient D can be estimated as the product of the cosmic ray velocity u ~ c and 
its mean free path Iq. A more precise analysis gives D = Iqv/S, where the factor three reflects 
the number of spatial dimensions. 



Dif Fusion coefficient We estimate now the energy-dependence of D(E) and its numerical 
value for a cosmic ray propagating in the Galactic disc. We start picturing its propagation as 
a random- walk caused by scatterings on magnetic clouds of size tq. Then one can distinguish 
two different regimes: 

• At low energies, i.e. when the Larmor radius Rl = p/{ZeB) is smaller than the size ro of 
magnetic clouds with density n, the angles between the entrance and the exit directions 
are isotropically distributed. Since the direction is on average changed considerably in 
each scattering process, the mean free path Iq is simply the distance between clouds, 
Iq = l/((jn) ~ l/{rQn) and thus 

1, 1 c , , 

-Co = -Iqv ~ — ~ const. (3.30) 
3 3 VqU 



• At high energies, cosmic rays are deflected in each cloud only by a small angle 5 ~ tq/Rl. 
The directions are uncorrelated and thus the mean deflection is zero, {6) = 0, and the 
variance (5^) is given again by the result for a random-walk, (5^) ~ N{rQ/ Rl)"^ . The 
effective free mean path Iq is the distance after which ((5^) ~ 1. Hence the energy- 
dependence of the diffusion coefficient is 

D{E) = {—y" Dq E"^ . (3.31) 



• The transition between these two regimes happens when RL(Ecr) = tq. Numerically, 
this energy is given by Ecr ~ 10^^eV(i?//xG)(ro/pc). 

Obviously, the picture of magnetic clouds or domains with an unique size ro is an oversim- 
plification. In a more realistic picture, there is a distribution of magnetic field fiuctuations 
that can be easiest characterized by the spectrum of its Fourier components, {B'^{k)) oc k^°'. 
Charged particles scatter mainly at field fluctuations which wave numbers k matches their 
Larmor radius, k ~ 1/Rl- If the amplitude of this resonant magnetic held fluctuation is (5i?res, 
then D {6Bres/B)-'^vRL/3 instead of our over-simplified estimate ()3.30p . Thus the energy 
dependence of D below Ecr is determined by the power-spectrum of magnetic field fiuctu- 
ations. The size ro of magnetic field domains is in this picture replaced by the correlation 
length Ic, i.e. the length scale below the field is smooth. 
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DifFusion and cosmic ray anisotropies In the diffusion picture, the resulting net cosmic ray 

flux is connected via Fick's law with the diffusion tensor. We assume a small, pure dipolc 
anisotropy, / = /q + /i cosi} and /i <C /q; and choose the z axis along the dipole axis. Only 
the dipole term contributes integrating I{E,'d), 

J^^{E) = 2n Msm^I{E,^)cos^ = —Ii{E). (3.32) 
Jo 3 

On the other hand, for diffusion Fick's law is valid, 

J^^{E) = -D^^dn{E)/dz . (3.33) 
Comparing the two equations we obtain 



A=-^Z>.|. (3.34, 



Then the anisotropy 6 is 



-'max + -'min lo 47r/o CU dz 

The experimental results for 5 provide tluis information on D. For an estimate we set 
dn/dz n/h where h is the characteristic scale on which n changes. With h ^ 500 pc and 
using ~ 3 X 10"*^, it follows 

D K\5ch^bx 10^^ cmVs . (3.36) 

Next we can estimate the mean free path of cosmic rays from D = cIq/3 and obtain Iq ~ 
5 X lO^^cm or 0.15 pc. Finally, we estimate the residence time of cosmic rays in the disc as 

r = ^«8xl0V- (3.37) 

This means that the prediction for the residence time in our simple diffusion picture agrees 
in magnitude with observations. Since cosmic rays are confined in the Galactic plane, and 
the Sun is located somewhat above the plane, this diffusion flux should be directed towards 
the northern Galactic hemisphere. 

Complete cascade equation A rather general set of equations describes the evolution of N 
species coupled by interactions k i are the following transport or cascade equations, 

^Ih^^^-^{DVni{E,^,t)) = Q{E,^,t) 

- {cpXrL{E) + X^')n^{E,^,t) 

- —iPm{E,:K,t)) 
^, dakiiE',E) 



+ X^J^dE — nk{E , X, t) . (3.38) 



k 

The first line describes diffusion, while the second line accounts for the loss of particles in the 
energy interval [E : E + dE] because of interactions or decays with = Tti/2- The third 
line describes continuous energy losses /? = dE/dt of a particle i: Two important examples 
are synchrotron radiation and adiabatic redshift losses due to the expansion of the Universe. 
Since f3n is the particle flux in energy space (with "velocity" dE/dt), d/dE{Pini{E,yi) is the 
divergence of this flux. Hence it is the analogue in energy space to V^j in coordinate space. 
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Leaky box model This model assumes that cosmic rays inside a confinement volume (e.g. 
the disc) have a constant escape probability per time, t^sc S> c/h. Neglecting all other effects, 

\. ' = = DAmiE, x) , (3.39) 

Cl- ''"esc 

and the distribution of cosmic ray path lengths in the disc is rii = no exp(— t/resc) = 
no exp(— z/Aesc)- Hence one can replace 

DAmiE, x) (3.40) 

in this model. Physically, the diffusion coefficient is a function of the distance to the disc, 
D = D(z), and the escape probability increases and the cosmic ray density decreases for 
increasing \z\, while in the leaky box model both are constant inside the confinement volume. 
Therefore, the leaky box model can describe only average values of observables that depend in 
general on z and care has to be taken comparing quantities calculated within the two different 
models. 

If we consider now the steady-state solution, drii/dt = 0, and replace the diffusion term by 
rii/Tesc, then 

For primary types like protons or iron, the decay term vanishes and production via fragmen- 
tation can be neglected. Introducing Aesc = PcpT^sc as the amount of matter traversed by a 
particle with velocity /3c before escaping, we obtain 

n- ~ Qi'^esc ^2 ^2) 

1 + -^esc/ \ 

The escape time Tesc in the leaky-box model should be, similar to D in the diffusion model, 
energy dependent. For the simplest hypothesis that Tesc{E) of different elements depends only 
on Z, one obtains from a fit to data 

g /4ZGeV\'' 

Aesc «i llAr for p > 4ZGeV (3.43) 

cm^ VP/ 

with 6 0.6 and Aesc = const, at lower energies. 

For protons, Ap = 55g/cm^ 2> Aesc for all energies, and thus 

rip = QpTcsc oc QpE~^ . (3.44) 

Hence the generation spectrum of protons should be steeper than the one observed, Qp oc 
^-2.7+d ^ ^-2.1 

For the other extreme case, iron, the interaction length is Ape = 2.6g/cm^. Hence at 
low energies, iron nuclei are destroyed by interactions before they escape, Ape *C Aesc, and 
therefore the iron spectrum reflects the generation spectrum, npe oc Qpe- Starting from 
the energy where Ape ~ Aesc, the iron spectrum should become steeper. The observed iron 
spectrTim is indeed flatter at low energies and steepens in the TeV range. 

A main test for propagation models are radioactive isotopes with life-time Ti/2 ~ "^esc- The 
abundance of such isotopes, e.g. ^''Be with ti/2 = 3.9 x 10^ yr, can be used to deduce separately 
Tesc and the density of gas. One finds as mean density of the traversed gas n ~ 0.3/cm^, i.e. 
just one third of the value in the disc, supporting the idea that cosmic rays are confined not 
only within the disc but in an extended halo. 
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Figure 3.7: The cosmic ray intensity in the knee region [9]. 

Knee The knee marks a break in the cosmic ray spectrum at ~ 3 x 10^^ eV. There are 
three classes of explanations: (i) a break in the diffusion coefficient as function of energy, 
(a) a signature of the -Emax distribution of sources, {in) a change in the interactions of 
cosmic rays in the atmosphere at ^/s ~few TeV. The possibilities (i) and (ii) produce both a 
rigidity-dependent knee, i.e. the position of the knee for different nuclei should be connected 
by Ez = ZEp, while for (in) the position of the knee depends on A. 

Experimentally, a rigidity-dependent knee is favored as it is visible from Fig. 13. 8[ The 
energy of the knee agrees roughly with the critical energy, when the Larmor radius of a 
proton becomes equal to the maximal length of magnetic field domains. Thus the knee might 
correspond to a transition from D oa E to D (x E^ , more generally to an increased leakage 
of cosmic rays out of the galaxy. 

3.3 Exercises 

1. Derive the solution of Eqs. (j3.22p and that the observed ratio 0.25 is reproduced for 
X PS 4.3g/cm^ 

2. Show that the diffusion equation with a point source at r = and a space-independent 
diffusion coefficient has the solution n{E,r) = l/{4TTr)[Q{E)/D{E)]. 
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Figure 3.8: The chemical composition measured by the KASCADE experiment in the knee 
region. Additionally, the position of a rigidity-dependent knee is for p, He, C and 
Fe indicated [lOj. 
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4.1 Sources of high energy cosmic rays 

The changes in the slope of the energy spectrum and in the composition from protons to iron 
in the knee region suggest that there are at least two different kinds of cosmic ray sources. It 
is natural to associate the bulk of cosmic rays with Galactic sources, while the highest energy 
cosmic rays have an extragalactic origin, produced by more powerful galaxies than ours. A 
further argument favoring an extragalactic origin of cosmic rays with the highest energies is 
the observed isotropy of their arrival directions on large scales that suggests a cosmological 
distribution of their sources. 

4.1.1 General arguments 

Energy argument favoring supernova remnants (SNR) The luminosity Lqr of Galactic 
cosmic ray sources has to fit the observed energy density pcR ~ 1 eV/cm^ of cosmic rays, 
taking into account their residence time Tgsc ~ 6 x 10^ yr in the Galactic disk. With Vd = 
nR'^h ~ 4 X 10^^ cm^ for R = 15kpc and h = 200 pc as volume of the Galactic disc, the 
required luminosity is Lcr = VdPcr/tcsc ~ 5 x lO^'' erg/s. 

In a successful core-collapse supernova (SN) around IOMq are ejected with w ~ 5x10^ cm/s. 
Assuming 1/(30 yr) as SN rate in the Milky Way, the average output in kinetic energy of 
Galactic SNe is LsN,kin ~ 3 x 10^^ erg/s. Hence, if the remnants of SNe can accelerate particle 
with efficiency 0(0.01), they could explain all galactic cosmic rays as it was suggested first 
by Ginzburg and Syrovatskii in the early 1960s. 

Hillas argument The Larmor orbit Rl = E/{ZeB) of accelerated particles has to fit inside 
the accelerator of size Rg, Rl = E/{ZeB) < Rg. For known magnetic fields and source sizes, 
one can constrain thus the maximal achievable energy as E'max = TZeBRg, as it is done in 
Fig. 14.11 for a compilation of potential cosmic ray sources. The Lorentz factor F introduced 
in -Emax accounts for a possible relativistic bulk motion of the source and is probably only 
for gamma ray bursts a significant correction. In Fig. 14. 1^ a so-called "Hillas plot," sources 
able to accelerate protons to E > 10^^ eV should lie above the solid red line, while sources 
above the green line can accelerate iron up to 10^'^ eV. Few candidate sources for acceleration 
to £^ = 10^1 eV seem to be compatible with Hillas' argument. 

On one hand, this constraint looks like a solid upper limit, because energy losses are ne- 
glected, the maximal acceleration time is finite and no accelerator is 100% efficient. On the 
other hand, nonlinear processes may lead to an amplification of magnetic fields inside the 
source. In general, sources neither too small (minimizing energy losses) nor too big (avoiding 
too long acceleration times) are favored. 
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Figure 4.1: Magnetic field strength versus size of various suggested cosmic ray sources. 

Blandford argument Tlie acceleration of a proton to the energy E = 10^^ eV by regular 
electromagnetic fields requires the potential difference U = lO'^^ V. What is the minimal power 
P dissipated by such an accelerator? In order to use the basic equation P = UI = U'^/R 
known from high-school physics, we have to know the appropriate value of the resistance R. 
Since the acceleration region is in most cases nearly empty, we use R ~ 1000 O (lead by the 
"impedance of the vacuum", R = Anko/c = l/(eoc) ^ 377 Q). Hence a source able to produce 
protons with E = W^^ eV by regular acceleration in electromagnetic fields has the minimal 
luminosity [12] 

L = C/V-R > lO^'^W = lO^'^erg/s. (4.1) 

This can be transformed into an upper limit on the density Ug of ultrahigh energy cosmic 
rays (UHECR) sources, since the observed UHECR intensity fixes the required emissivity £, 
i.e. the energy input per volume and time, as £ ~ 3 x 10^^erg/(Mpc^yr). Hence, the density 
of UHECR sources able to accelerate protons to E = W^^ eV should be smaller than = 
C/L 10~^/Mpc^, if the acceleration is by regular electromagnetic fields. For comparison, the 
density of normal galaxies is Ug ~ 10~^/Mpc^, while the most common type of active galactic 
nuclei in the nearby Universe, Seyfert galaxies, has the density ~ (1 — 5) x 10~^/Mpc^ 
within redshift z < 0.02. 

4.1.2 Specific sources 

Most galactic astrophysical sources are connected with type 11 (or core-collapse) supernovae 
(SN) and their remnants (SNR): Examples are the direct acceleration in the magnetosphere of 
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young pulsars and shock acceleration in SNRs. Above E > 10^^ eV, the signature of galactic 

sources would be an enhanced flux towards the galactic plane. Practically all extragalactic 
sources except gamma ray bursts are associated with active galactic nuclei (AGNs). 



SNe II Type II or core collapse supernovae occur at the end of the fusion process in very 
massive stars, M > {5—8)Mq. Theses stars develop an onion-like structure with a degenerate 
Fe core. After the core is completely fused to iron, no further processes releasing energy are 
possible. Instead, photo-disintegration destroys the heavy nuclei, e.g. via 7+^^Fe —>■ ^He-|-4n, 
and removes the thermal energy necessary to provide pressure support. In the following 
collapse of the star, the density increases and the free electrons are forced together with 
protons to form neutrons via inverse beta decay, e~ + p ^ n + Uc.: A proto- neutron star 
forms. When the core density reaches nuclear density, the equation of state of nuclear matter 
stiffens and infalling material is "reflected," a shock wave propagates outwards heated by 
neutrino emission from the proto-neutron star. If the SN is successful, a neutron star is left 
over; otherwise a black hole remains. 

The released gravitational binding energy, 

is emitted mainly via neutrinos (99%). Only 1% is transferred into kinetic energy of the 
exploding star and only 0.01% goes into photons. 
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Pulsars are the left-over of the former iron core, now transformed into a neutron star with 
Rns/Rq ~ 10~^. They may be born fast rotating with a strong magnetic field, because of 
conservation of 

• angular momentum, L = lu) = const., with / ~ MR?: A star like the sun, rotating 
once per month, P 10^ s, would rotate with a ms period when contracted down to 10 
km in size. 

• magnetic flux, (pB = BA = const., for an ideal conductor, where magnetic field lines 
are frozen in the plasma. As the core collapses, the magnetic field lines are pulled 

more closely together, intensifying the magnetic field by a factor (i?0/i?Ns)^ ^ lO^''. 
Magnetic A stars have surface fields up to 10^ G, while the maximal observed field 
strengths of neutron stars are of the order B ~ 10^^ G. 

Rotating dipole model The energy of a rotating sphere \s E = luP' /2 and the energy change 
is 

E = IujCo . (4.3) 

The rotational kinetic energy of the Crab nebula, the remnant of a SN observed by Chinese 
astronomers in 1054, is E = \Iuj^ « 3 x 10^^ erg with I = aMR^ and a « 2/5 for a 
homogeneous sphere; its energy-loss per time is the time-derivative E = luu ^ 7x 10^^ erg/s, 
with 7 = 1.5 X 10^^ gcm^ and ci; = 4 x 10~^yra;. 

If the energy is lost as electromagnetic dipole radiation, then 

w C-"' 
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where a denotes the angle between rotation and dipole axis, and we obtain an estimate for 
Crab's magnetic field, i? ~ 7x 10^^ G for sin a = 1. Thus pulsars have indeed extremely strong 
magnetic fields that, if they are fast rotating, may accelerate particles to high energies. 



Direct acceleration of particles by pulsars The light-cylinder around a pulsar is the surface 
at Rc = c/uj. If the magnetic field lines would rotate rigidly with the pulsar, then the linear 
velocity reaches the speed of light, v = c, at Rc- 

If the magnetosphere is filled sufficiently with plasma, the electric conductivity o" ^ oo, 
and Ohm's law J = a{E + vB/c) implies 

E = -vxB/c, (4.5) 

where v = a; x r. 

For a magnetic dipole field pointing along = 0, 

B(r) = {2er cos + sin t?) , (4.6) 

the potential difference between the polar cap and infinity follows as 

A(/> = y E • ds = y V X B • rfs = ^"^^'^ sin^ i9 . (4.7) 

Open field lines that extend beyond the light cylinder c/iv start near the polar cap, within 
the cone ??o ~ \/ioR/c. Inserting i?o into Eq. (j4.7p we obtain as maximal acceleration energy 
achievable by a pulsar 

c c2 1013 G V 3000 s-i/ ^ ^ 

Thus a young, fast rotating pulsar appears to be a very good particle accelerator. The 
main problems are that in realistic models the potential difference Acj) that can be used for 
particle acceleration is much smaller and the extreme energy losses due to, e.g. curvature 
radiation. The magnetosphere of a pulsar may be also filled by 7 ^ e~^e~ with a plasma. 
Finally, pulsars as main source of UHECRs would predict a strong anisotropy of the UHECR 
intensity, because neutron stars are concentrated in the Galactic plane. 



Shock acceleration in SN remnants is the standard paradigm for cosmic ray acceleration 
up to ~ lOi^ — lOi'^eV in our galaxy and is discussed in the next section. 



Active galaxies is a common name for all galaxies with unusual emission that is not asso- 
ciated with stars. In contrast to normal galaxies whose total luminosity is the sum of the 
thermal emission from each of the stars found in the galaxy, a large fraction of the total 
luminosity of an active galaxy is non-thermal and is emitted by the nuclei of the galaxy. 
The modern view is that the common mechanism behind the energy generation in AGNs is 
accretion on the SMBH in their center. 

Main reason for this view is the fast variability of their spectra, from time-scales of months 
for quasar spectra down to daily variations for blazars. Causality limits then the size of 
the emission region as ct < 200 AU. Secondly, the enormous energy output of AGNs requires 
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Figure 4.2: The unified scheme of AGN. 

an extremely efficient energy generation mechanism. For accretion on a BH, the maximal 
energy gain is -Emax ~ GmM/Rs, where the Schwarzschild radius is Rs = 2GM/c^ , and thus 
-Emax = m(? /2. A large part of this energy will be lost in the BH, while the remainder heats 
up via friction an accretion disc around the black hole. Modeling the accretion process gives 
an efficiency of e = 10%-20%. Thus the luminosity from accretion is 

L='-^. (4.9) 

For a rather modest mass consumption of the BH, dm/dt = IMg/yr, one obtains 6x lO^^erg/s 
or L ~ 10^2^0. 

The unified picture of AGNs is illustrated in Fig. 14.21 The different AGN types are only 
facets of the same phenomenon-accretion on a SMBH — viewed from different angles, at differ- 
ent stages of activity (small or large dm/dt) and evolution in time (e.g. from a quasar phase 
at redshift z ~ 1-4 towards a Seyfert galaxy at present). Blazars are AGN with a relativistic 
jet that is pointing in the direction of the Earth and are therefore often ranked among the 
most promising sources of ultrahigh energy cosmic rays. According Fig. blazars are either 
Flat Spectrum Radio Quasars (FSRQ) or BL Lac objects. 

Gamma Ray Bursts come in two (main) sub- varieties, depending on their duration that 
varies from fraction of a second to many minutes. While short duration GRBs are most likely 
the result of binary collisions between e.g. neutron stars, long duration GRBs which make 
up about 2/3 of all GRBs are associated with supernova events in extremely massive stars. 
GRBs are highly beamed sources of gamma-rays and perhaps also of high energy neutrinos 
and cosmic rays. A distinctive feature is the high Lorentz factor of shocks in GRBs. 
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Figure 4.3: A cosmic ray "scattering" elastically at a magnetic cloud moving with velocity 
V. 



4.2 Acceleration of cosmic rays 
4.2.1 Second order Fermi acceleration 

We consider a cosmic ray with initial energy Ei "scattering" elastically on a magnetic cloud 
that moves with velocity V <^ c. We want to derive the energy gain = (E2 — Ei)/Ei per 
scattering. The variables and {E2,'&2) are shown in Fig. 14.31 which we label with a 

prime in the cloud system and without a prime in the lab system. 

With a first Lorentz transformation between the lab and cloud system we connect the 
variables characterizing the cosmic ray entering the cloud, 

E[=-fEi{l- Pcos^i) where (3 = V/c and 7 = 1/^1^^, (4.10) 
and with a second Lorentz transformation the exit variables, 

= 7^^2(1 + /3 cos t?^. (4.11) 



Since the scattering off magnetic irregularities is collisionless and the cloud is very massive, 
energy is conserved, E2 = E[. Hence we can eliminate E'2 and obtain as relative energy gain 

E2-E1 1 - /? cos 1?l + /? cos ??2 - cos 1?l cos T?2 n 

^ = — = Trp2 1- (4-12) 

To proceed, we need average values of cosi9i and cos??2- Since the cosmic ray scatters off 
magnetic irregularities many times in the cloud, its exit direction is randomized, (cos 192) = 0. 
The collision rate of the cosmic ray with the cloud is proportional to their relative velocity 
(v — ycos"!?!). For ultrarelativistic particles, f — > c, and the collision rate is 

n?7 

(X (l-/?cosi?i). (4.13) 



32 



4.2 Acceleration of cosmic rays 



We obtain (cosi?i) averaging cost?i weighted by dn/dJ7i over all angles, 

(,os<),>=/cos<),i^d«,//i^dn, = -f. (4.14) 

Plugging (cos ■i?2) = and (cos i^i) = — /3/3 into Eq. (|4.12|) and taking into account that /? ^ 1 
gives as average gain 

Thus {(,) oc > and we have shown that on average a cosmic ray gains energy scattering 
on "magnetic clouds" with an ordered bulk velocity V. The energy gain per scattering is 
however only of second order in the small parameter /3, and acceleration is therefore rather 
inefficient. Moreover, one can show that the resulting energy spectrum would depend strongly 
on the cloud parameters, making thereby the observed feature-less power spectrum of cos- 
mic rays difficult to understand. Acceleration at shocks that we consider next avoids both 
disadvantages. 

4.2.2 Shock or first order Fermi acceleration 

Ideal fluids and shocks As first step in gaining some understanding of shocks we recall the 
the ideal fluid equations. The three basic equations are the conservation equation for mass, 
for momentum, and the Poisson equation, 

dtp + V ■ (pv) = 0, (4.16a) 

P^ + Z'v- (Vv) =F- VP, (4.16b) 
at 

A(/. = 4^G/9. (4.16c) 

The LHS of (I4.16bl) . the Euler equation, measures the velocity change dv/dt of a fluid element, 
summing up the change at a fixed coordinate, dtv, and the change due to the movement of the 
element, while the RHS consists of an external force F and a force due to a pressure gradient 
VP. The Poisson equation connects the mass density p with the gravitational potential (p. 

We consider now small perturbations xi around a static background, po = const., Pq = 
const, and vq = 0. Restricting us to small perturbations, xi <C xq, allows us to neglect 
all quantities quadratic in the perturbations. Sound waves propagate in most circumstances 
adiabatically, i.e. without production of entropy, dS = 0. Thus changes in P and p are 
connected via 



dP\ fdP 



Inserting x = xq A- xi into the fluid equations and neglecting quadratic term gives 

+ poV • VI = , (4.18a) 

A01 = 47rGpi , (4.18b) 

dt^i + ^Vpi + V01 = , (4.18c) 
Po 
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where we used also dxP = Pi = -§^(?s- These three equations can be combined into 

one second-order differential equation for pi. We multiply Eq. (|4.18cp by pQ and apply V on 
it, 

cl^Pi = -po{dtV -^i + Mi)- (4.19) 

Then we insert ()4.18bp for /S.(f)i and (j4.18ap for V • vi, and obtain a linear, inhomogeneous 
wave equation 

d^pi - c^Api = 47rGpipo (4.20) 

pressure grav. force 

for the density perturbation pi. The dispersion relation of plane waves exp{—i{ujt — kx)), 

confirms that Cg = {dP / dpi)^^"^ is the sound speed. 

For a mono-atomic gas, the equation of state is P = Kp^ with 7 = 5/3. Thus, the sound 
speed is Cs = [jP/ p)^/"^ and if an adiabatic compression with density p2 = epi propagates, 
then Cs oc e^'''^^^/^ Hence the sound speed increases for a compression, the dense region over- 
runs uncompressed regions and becomes even denser: A discontinuity (=" shock") develops 
in some hydro dynamical variables like the density. 



Mach number of a (strong) shock Our main aim is to derive the Mach number M. = v/cg 
of a shock that as we will see determines the slope of the energy spectrum of accelerated 
particles. We will consider only the properties of an one-dimensional, steady shock in its rest 
frame and assume that magnetic or gravitational fields can be neglected. Then the continuity 
equation for mass, dtp + V • (pv) = 0, becomes simply 



da; 



(H =0. 



(4.22) 



The Euler equation simplifies using the same assumptions and taking into account Eq. (j4.22p 
to 

^{P + pv"^) = . (4.23) 



Additionally to the conservation laws for mass (|4.22p and momentum (|4.23p we need the 
conservation law for energy. 



P ' 



(4.24) 



Here, the first bracket accounts for the change of kinetic, internal and potential energy with 
time which has to be balanced by the energy flux through the boundary of the considered 
volume. Specializing again to the case of an one-dimensional, stationary flow with <1? = 
gives 




+ {U + P)v \ = Q. 



(4.25) 
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Figure 4.4: Conditions on the down-stream (left) and the up-stream (right) side of a shock 
in the lab system (top) and in the shock rest frame with vi = — (bottom). 



Integrating these equations over the discontinuity of the shock results in the "Rankine- 
Hugoniot" jump conditions, 



pv 



P + pv 

7 



Pv 



0, 
0, 

^0, 



(4.26a) 
(4.26b) 

(4.26c) 



where we used also [/ = P/(7 — 1). Since we assumed a steady flow, these jump conditions 
have to be evaluated in the shock rest frame, cf. Fig 14.41 - otherwise time-derivatives should 
be included. 

Inserting first p2 = {vi/v2)pi into (|4.26b|l gives P2 = Pi + pivi{vi — V2)- Next we use these 
two expressions to eliminate p2 and P2 from Eq. ()4.26c|) . Reordering the resulting equation 
according to powers of V2, 



7+1 
7-1 



2^ 27 
V2 + 



7-1 



^1 + PlVl 
pivi 



V2+V 



27 Pi 
7 - 1 pi 



0. 



(4.27) 



replacing Pi by the sound speed and dividing by vf, we obtain finally a quadratic equation 
for t = V2/V1, 



7 + 1 

7-1 



2 27 

H — 

7-1 



3+7 ^+1 + 



-f-lvf 



0. 



(4.28) 



Now we recognize vi/ci as the Mach number M. Since we are interested in fast flows, vi S> ci, 
we can neglect the two terms and obtain as approximate solutions 



t = 1 



or 



Vl =V2: 



(4.29a) 
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t = = R or Rv2 = vi . (4.29b) 

7+1 

The first solution is obviously trivial, while the second one is the strong shock solution. The 
compression ratio R indicates how strong the density and the velocity in the up- and down 
stream regions differ. Since we are using as reference frame the shock frame, we have Vs = vi 
and thus for 7 = 5/3 and R = 4, 

V2 = Vs/R = Vs/4, (4.30a) 
P2 = Rpi = 4/91 , (4.30b) 
P2 = ^pivl/4. (4.30c) 

Hence no matter how strong a shock is, it can compress a mono- atomic gas only by a factor 
four. The same factor relates the velocity of the shock and of the matter after the passage 
of the shock. It is this universal ratio of the up- and down-stream velocities for any strong 
shock that results in the generic prediction of an 1/ E"^ spectrum produced by acceleration at 
shocks. 



Acceleration at shocks We consider again a cosmic ray with initial energy E\ "scattering" 
at magnetic irregularities that are now separated by a planar shock moving with velocity Vs- 
The only change to the previous discussion are the resulting different angular averages. The 
(normalized) crossing rate is given by the projection of an isotropic flux on the planar shock, 

dn _ f 2cos??i cos'!?i<0 ^Al.^\ 
dcos??i ~ \ cosi?i > ■ ^ ^ 

while the crossing rate dn/dcosi?2 is non-zero for cost?2 > 0. Thus (cost?i) = — | and 
(cos'!?2) = I and hence 

The average gain {^) is now linear in (3 and justifies the name "first-order" Fermi (or shock) 
acceleration. 



Energy spectrum produced by Fermi acceleration The energy En of a cosmic ray after n 
acceleration cycles is 

£^„ = ^0(1 + 0" (4.33) 
and the number of cycles needed to reach En is thus 

n = ln(^)/ln(l + 0. (4.34) 

If the escape probability Pesc per encounter is constant, then the probability to stay in the 
acceleration region after n encounters is (1 — Pesc)"- We obtain the fraction / of particles with 
energy > En as 

/(> E)=± (1 - Peso)- = ^^^^^ oc -L , (4.35) 

m='n Peso Peso \^0/ 



36 



4.2 Acceleration of cosmic rays 



where 

7 = In (t^—) I + « V.Ji (4.36) 

VI -PescJ 

with ^ ^ 1 and pesc ^ 1- Hence, both first and second order Fermi acceleration produces a 
power-law energy spectrum. 

The maximal energy achievable by a specific source is determined by several factors. First, 
the finite life-time limits the number of cycles n and thus Second, the escape probability 
is energy dependent and increases generally for increasing energy. Third, energy losses like 
synchrotron radiation increase with energy and balance at a certain point the energy gain. 

Exponent 7 for shock acceleration Since we know already that ^ f« |/9 = | ^^~^^ , we have 
to determine only the escape probability pesc in order to estimate the exponent 7 Pesc/i 
of the integral spectrum produced by shock acceleration. The particle flux T through an 
infinite, planar shock front is (cf. Eq. (j3.2p ) 

HE) = TTliE) = ^ , (4.37) 

assuming Vg <^ c and an efficient isotropization of cosmic rays up-stream. 

In the shock rest frame, there is a particle flow ^esc(-E') = V2n{E) downstream away from 
the shock front that will be lost for the acceleration process. Thus the escape probability pesc 
as the ratio of the loss and the crossing flux is 

J^esc V2n 4V2 , . 

Pesc = —p- = — 77 = (4.38) 

J- cn/4 c 

and the spectral index of the integral energy spectrum follows as 

7 ~ Pcsc/e ~ -j^—r . (4.39) 

V1/V2 - 1 

Typical values for the sound speed in the interstellar medium around a SNR are ci ~ 10 km/s, 
while the shock velocities are vi ^ 10^ km/s. Thus the Mach number Ai = vi/ci ^ 1 and 
we can use the result vi = Av2 in the strong shock limit. The exponent predicted by Fermi 
acceleration at non-relativistic shock is therefore independent of the shock parameters and 
agrees with the value needed to explain the spectrum of Galactic cosmic rays. These are the 
two main reasons for the popularity of shock acceleration. 

We can check that possible corrections are small by solving Eq. (j4.28p for arbitrary M. 
Then we obtain R = AM'^/{'i + M'^) and thus 

Clearly, other effects like energy losses or an energy dependence of pesc will have a larger 
impact on the resulting spectrum than the 0{M.~'^) corrections. Moreover, we worked within 
an extremely simplified picture and should expect therefore deviations from a power-law with 
7 ~ 1. For instance, we used the test particle approach that neglects the influence of cosmic 
rays on the plasma. 

Finally we remark that the observed flux from extragalactic sources might require a steeper 
injection spectrum than 7 ~ 1, and that the acceleration mechanism responsible for cosmic 
rays beyond ~ 10^^ eV is still open. A discussion of more recent developments in the under- 
standing of cosmic ray acceleration can be found in Ref. |13j . 
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5.1 Experiments and detection methods 

The Earth's atmosphere is opaque to photons with energy above 10 eV, meaning that to 
observe gamma rays directly requires placement of a detector above the earth atmosphere. A 
major turning point in gamma-ray astronomy was therefore the launch of the first satellite- 
borne telescope, SAS-2, in 1972. 

The heretofore most successful gamma-ray satellite was the Compton Gamma Ray Observa- 
tory, taking data from 1991-2000 with its four different experiments. Some of the main results 
were the proof that a large fraction of gamma ray bursts are isotropically distributed and thus 
extragalactic in origin, the detection of discrete sources of extragalactic 7-ray emission, and 
an upper limit on the diffuse extragalactic 7-ray background. A new satellite experiment, 
GLAST, will be launched 2008. 

Main limitation of satellite experiments is their small collection area, which limits their use 
to energies < 100 GeV. On the other hand, 100 GeV is the limit where the electromagnetic 
cascade in the earth's atmosphere from the initial photon can be detected. While the cascade 
dies out high in the atmosphere below 10 TeV, showers are still detectable via the Cherenkov 
emission of relativistic electrons and muons. The main difficulty in ground-based gamma-ray 
astronomy is the presence of an almost overwhelming background of charged cosmic rays. 

Discrimination The image of an air shower in an atmospheric Cherenkov telescope can be 
modeled as a two-dimensional ellipse. The two main parameters that distinguish an air 
shower initiated by a photon and a proton are the length and the width of this ellipse: While 
a photon shower is narrow with a long elongation, a proton shower of comparable energy 
is much wider and shorter, cf. Fig. 15.11 The resulting discrimination power would be not 
sufficient to separate photons from the cosmic ray background, if photon sources like e.g. 
blazars would not appear point-like. 

5.2 Electromagnetic processes 
5.2.1 Synchrotron radiation 

The power P of emitted synchrotron radiation by a particle with mass m, charge e and 
momentum p± perpendicular to the magnetic field is in the classical limit 

^^' = 3""^ [m^) =3""^ [mBTj ■ ^'-'^ 

The energy loss per time is (3 = dE/dt = —P oc and hence is most severe for the 

lightest charged particle, the electron. The typical length-scale Zgyn of synchrotron losses 
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Figure 5.1: The longitudinal development (top) and the projection (bottom) of a shower 
initiated by a photon (left) and by a proton (right) [11]. 
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{Bcr = rrig/e = 4.41 x 10^^ G is the critical field strength for an electron) 

_ / I dEy^ _ 3 1 fB,A' 

For electrons with E = 10^^ eV in the Galactic disc, B ~ 3/xG, we have /gyn ~ 100 pc. Thus 
one of the reasons why we have not discussed electrons as a major component of cosmic 
rays is that their energy losses are too severe. On the other side, synchrotron radiation of 
relativistic electrons is in the radio range and explains the strong radio emission of many 
AGNs: Synchrotron radiation peaks at 0.29wc with 

^c = i7^rn. (5.3) 

For B ~ 0.1 G close to a source and 7 = 10^, the peak of synchrotron radiation is at 
lO^^'^me ~ 100 GHz. Synchrotron radiation is the explanation of the first peak visible in 
Fig. 15.21 showing the observed spectrum from a Blazar. 

5.2.2 (Inverse) Compton scattering 

The cross section for Compton scattering e~ + 7 — > e~ + 7 is with crxh = Svra^/ (3mg) 

o" = 0-Th (1 - s/m^ + . . for s/m^ <C 1 (5.4) 

3m? A / 2 1 An / 2 

a = — — CTTh Ins/m H \- . . . ] for s/m »1 (5-5) 

4s V 2 / 

i.e. Compton scattering is suppressed for s » m^. 

While in the classical Compton experiment an energetic photon hits an electron at rest and 
transfers part of its energy, in astrophysics "Inverse Compton scattering" is more important: 
A fast electron hitting a low-energy photon transfers a large fraction y of its energy to the 
photon, 

y w l/ln(s/m^) . (5.6) 

This is the explanation of the second peak visible in Fig. 15.21 showing the observed spectrum 
from a blazar. 



5.3 Hadronic processes 

Neutral pion decay vr'^ 27 In the rest system of the vr'^, the two photons are emitted 
back-to-back, pi = — p2 and Ei = E2 = m-,j/2. How are the photon energies distributed, if 
the pion was moving with velocity vl 

Let us use the Lorentz transformation between energies in the rest and the lab system. In 
the lab system, 

E = -i{E' + (3p' cos 'd) (5.7) 

where •& is the angle between the velocity [5 of the pion and the emitted photon. The maxi- 
mal/minimal values of E follow directly as E = ^{E' it [5p') for cos = ±1, i.e. if the photons 
are emitted parallel and anti-parallel to the direction of flight of the pion. 
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Inserting E' = p' = gives 

E^n = \ {irn^ ± (iim.) = \E^{l±f3) (5.8) 

or -Emax -E-TT and Smin — > in the ultra-relativistic limit /?—>!. 

What is the distribution of the emitted photons? Since the pion is a scalar particle, the 
photon distribution is isotropic in the lab system, 

diV = ^dSl = ^d| cos d\ . (5.9) 

Using (j5.7p . one can express dE as 7/3p'd| cosi?| and thus 

1 (\E 

AN = —dVt = . (5.10) 

Hence dA^/d-E = const., and the photon spectrum produced by a pion beam with uniform 
velocity is a box between -Bmin and -Emax from Eq. (j5.8p . If we use log(S) as x coordinate, 
the boxes are symmetric with respect to 7X1^^12. 



Photon spectrum from " many pion decays" Consider now the photon spectrum not from a 
single pion, but from a beam of pions with arbitrary energy spectrum. Each photon spectrum 
from a single decay is as we saw a box centered at m7r/2. Hence also the total photon spectrum 
from pions with an arbitrary energy distribution is symmetric with respect to 



5.4 Source models 

The energy spectrum observed from blazars, i.e. AGN which jets are pointing towards us, 
covers all the electromagnetic spectrum from radio frequencies up to TeV energies. It is char- 
acterized by two bumps, cf. Fig. 15.21 one peaking between the IR and the X-ray band, another 
one at gamma-ray energies. The first peak can be explained as synchrotron radiation from 
relativistic electron, while the second one is presumably due to inverse Compton scattering - 
between the same electrons and either soft photons or synchrotron photons. 

Although sources like blazars should accelerate not only electrons but also protons, syn- 
chrotron radiation and inverse Compton scattering of electrons are sufficient to explain in 
most cases the observed spectra. There are only few exceptions where an additional photon 
component from pion decay at the highest energies might be needed. 

The proof that a source accelerates indeed protons might be easier for Galactic sources, 
where the source region can be resolved by modern atmospheric Cherenkov telescopes. An 
example is the recent observation of the SNR RX J1713. 7-3946 by the HESS array that 
has been interpreted as the first direct evidence that SNRs accelerate cosmic rays up to 
TeVs |14j . The image of this source shows an increase of TeV photons in the directions of 
known molecular clouds. This is in line with the expectation that, if photons are produced 
as secondaries in pp collisions, their flux is correlated with the matter density in the SNR. 
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Figure 5.2: Spectrum from the blazar Markarian 421 in a flaring and a quiet state [15 



5.5 Electromagnetic cascades 

A high energy particle interacts at the top of the atmosphere and initiates a cascade. In 
the 1930's, it was recognized that the observed cosmic rays on ground are only the secon- 
daries produced in this cascade. Heitler introduced the following simple toy model of an 
electromagnetic cascade: 

• In each interaction characterized by the constant interaction length A, a parent particle 
splits into two daughter particles. Each particle takes half of the energy. 

• After n = X/X generations, the cascade consists of N{X) = 2^1^ particles. 

• The energy per particle is E{X) = Eq/N, where £"0 is the energy of the primary. 

• Particle production stops, when E{X) < E^r ~ 0.1 GeV. Below this energy, photons 
loose energy mainly by ionization losses. Unstable particles decay and the number of 
particles in the shower decreases. 



The number of particles at the shower maximum X^ax is 



N{X,, 



Eq 

En 



with Xr, 
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In- 



E, 







In 2 Er 



Thus the energy Eq of the primary particle is connected via 



N{X„ 



oc Eo 



and Xrriax oc In Ei 







(5.11) 



(5.12) 



to measurable quantities. Note that Ecr — the quantity that mainly controls N{X„ 
determined by well-known low energy processes. 
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Cascade equations and scaling The cascade equation for electron (= e^,e ) and photons 



are 



- nai^tot{E)Ni{E,r) +n dE — Nk{E ,r) (5.13) 

Je oE 



dr 

where i = 6,7. Two possible boundary conditions are 

iV,(^,0) = 5{E-Eq), (5.14a) 
Ni{E,^) = KE-", (5.14b) 

i.e. the study of a single shower or the general properties of electromagnetic cascades. We 
consider next only one species (i.e. we are not interested which particles are produced in the 
interactions) and the simpler second case for the boundary condition. Under which conditions 
may we hope for an analytical solution? The energy and depth dependence should factorize, 
N{E,r) = g{E)f{r), i.e. the shower should have at each depth the same energy dependence. 
This happens, if (z) the cross sections are not energy dependent, and (ii) the differential cross 
sections scale, 

a dz a oE 

with z = E/E' . Then, we can rewrite N{E',r) = KE'~" = Kz°'E~°' for boundary condition 
(b) and it follows 

N{E, r) = N{E, 0) exp[-nar{l - Z{a)] (5.16) 

with the spectrally averaged energy loss fraction 

ri 
10 



Z{a) = / dzz'^-'fiz). (5.17) 



Equation (I5.16P is straightforward to interpret: For z — > 0, one obtains the limiting case of 
full absorption, while interactions without energy transfer do not change N{E, r) at all. In 
the intermediate case, Eq. (|5.17p weights the energy loss z properly not only with f{z) but 
also with the spectral shape N{E) oc E~°^. 

For a steeply falling spectra, a > 1, the region z ^ 1 gives the dominating contribution to 
the integral defining Z(a). This illustrates that for the modeling of air shower the knowledge 
of cross sections in the "forward" region z ^ 1 is crucial. 

5.5.1 Cascades in the atmosphere 

Probabilities The differential probability P{z) for an electron with energy E scattering on 
an atomic nucleus to emit a photon with energy zE traversing the depth dt = dX/Xo is in 
the high energy limit 

Pey{z)=Z + ^^(^^ + 2b^ , (5.18) 

where 6 = b{Z) ~ 0.012 depends on the average charge Ze of air. This process is infrared- 
divergent, i.e. the probability to emit a soft photon diverges, Pe-y{z) — > 00 for z — > 0. Only if 
we consider the energy loss rate, 

§ = -^/'ci..EP„,.) = -S^.-f . (5.19) 
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we obtain a finite result. By contrast, the probability for the pair production of a massive 
e~^e~ pair in the Coulomb field of a nucleus is infrared-finite, 



^7e(^) = ^-^&+(^ + 2fo) (^-^) , (5.20) 



and can be integrated to 



J- = |d.P,.W = IJ.I. (5.21) 

Both electrons and photons lose their energy exponentially, E{X) = E'q exp(— X/Xq), with 
the radiation length Xq 37.2g/cm^ (in air) and 7Xo/9 as scales, respectively. Below the 
critical energy E^r ~ 84.2 MeV, ionization losses become more important than pair production 
and the growth of the cascade stops. 

Shower profile in the atmosphere The solution of the coupled cascade equations for elec- 
trons and photons in the high-energy limit, i.e. using Eq. (j5.18p and (j5.20p as splitting 
probabilities, and for a power-law boundary condition ()5.14bp at the top of the atmosphere 
is lengthy but straightforward. Here we mention only the two essential steps, for details see 
the classic paper of Rossi and Greisen [16]: First, one obtains an infrared- finite expression 
writing the divergent total bremsstrahlung probability as an integral and combining it with 
the gain term. Second, in the scaling limit the functions Z are for a fixed a simple numbers 
and one has to solve therefore only two coupled first-order differential equations. Hence, the 
solutions are of the type Aexp(Ai(X)) + Bexp(A2(X)). 

The number of electrons A'^e(^) in an electromagnetic shower initiated by a photon as 
function of the shower age 

s = (5.22) 



t + 2f3 



is approximately 



.r 0-31 

iVe = -^exp 



t ( 1 - ^ In s 



(5.23) 



where t = X/Xq and /? = In^E/Ecr)- This approximation takes into account also low-energy 
processes; Ne{X) is shown for several shower energies in Fig. 15.31 A shower initiated by a 
primary with energy E = 10^^ eV requires about one atmospheric depth to develop fully and 
contains around 10^^ particles at its maximum. 

5.5.2 Cascades on diffuse photon backgrounds 

A high energy photon scattering on a low energy background photon (e.g. from the cosmic 
microwave or infra-red background) can start an electromagnetic cascade, via the two main 
processes 

7 + 76 ^ e+ + e" (5.24) 

e + 7b^e+7. (5.25) 

Pair production stops for s = 4i?^e-y < 4mg, where is the typical energy of the background 
photons. Electrons continue to scatter on cosmic microwave photons in the Thompson regime, 
producing photons with average energy 
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Figure 5.3: Number of electrons Ne{X) in an electromagnetic shower initiated by a photon 
as function of the depth X for primary energies 10^^, 10^^, . . . , lO^'^ eV; the shower 
age s is also indicated. 

Hence the Universe acts as a calorimeter for electromagnetic radiation, accumulating it in the 
MeV-GeV range. Moreover, it means that the Universe is opaque for high energy photons. 
Figure 15.41 show the absorption length of photons on various sources of background photons. 
Starting from energies in the TeV range, high-energy photon astronomy can study only the 
local (and therefore also recent) Universe. This is one of the main motivations to explore also 
high energy neutrinos, despite of the experimental challenges in neutrino detection. 

5.6 Exercises 

1. Show that Pel for synchrotron radiation has the correct Lorentz transformation proper- 
ties by expressing it as function of F^^, and p^. 

2. Rewrite the energy losses due to synchrotron radiation and Compton scattering as 
function of the energy density of the magnetic field and real photons and show that the 
two expressions agree in the classical limit. 

3. Find the number of particles in the shower maximum, the corresponding shower age s 
and depth Xmax from Eq. ()5.23p . 
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Figure 5.4: The gamma and the cosmic ray horizon, i.e. the free mean path of photons and 
cosmic rays as function of energy: The (red/blue) shaded area is inaccessible for 
photons/cosmic rays. The distance to the nearest galaxy cluster is also indicated. 
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6.1 Propagation of extragalactic particles 
6.1.1 Energy losses of protons and nuclei 

There are three main energy loss processes for protons propagating over cosmological dis- 
tances: Adiabatic energy losses due to the expansion of the universe, —{(iE/dt)/E = Hq, 
e^e~ pair-, and pion-production on photons of the cosmic microwave background (CMB). 

The relative energy loss per time of a particle (due to interactions with the CMB) can be 
estimated as 

1 , ^ / ^ 

^ ^ = {y^n,) (6.1) 

where y = {E — E') / E is the energy fraction lost per interaction, ~ 410/cm^ is the density 
of CMB photons with temperature T ~ 2.7 K and the brackets (...) remind us that we should 
perform an average of the differential cross section with the momentum distribution n(p) of 
photons. In our estimates, we avoid this complication considering only reactions well above 
the threshold energy when essentially all photons participate in the reaction. 



Electron-positron pair-production Since the produced e"*'e~ pair in the process p + 7 — > 
p + + e~ is light, this energy loss process has a low threshold energy but leads in turn only 
to a small energy loseS per interaction, y = 2me/mp ~ 10 ^. The threshold energy on CMB 
photons follows from {kp + k^)"^ > (mp -|- 2me)^ as Ep > nif.'mp/E^ ~ 2 x 10^^ eV. The cross 
section of the reaction is o" = a/32(Txh/('S), i.e. is a factor ~ 2 x 10~^ smaller than Compton 
scattering. 



Pion-production and the GZK effect Greisen and Zatsepin and Kuzmin noticed that pion- 
production should introduce a "cutoff" in the cosmic ray spectrum. Since y ~ ra-,^/rap ~ 0.2 
close to threshold, but the cross-section is of similar size as the one of pair-production, the 
energy losses of protons suddenly increase around 5 x 10^^ eV. 

For an estimate of the cross section for pion production we use the Breit-Wigner formula 
with the lowest lying nucleon resonance as intermediate state, 

- (2,si + l)(2., + l) ^ {E-Ma? + TVA ' ^^-'^ 

where p^^, = {s - m%f/{As), Mr = tua = 1.230 GeV, J = 3/2, Ttot = 0.118 GeV, = 
0.55%. At resonance, we obtain ctbw ~ 0.4mbarn in good agreement with experimental data, 
cf. Fig. 16. 1[ We estimate the energy loss length well above threshold i?th ~ 4 x 10^^ eV with 



^ Close to the threshold, the energy fraction lost can be obtained considering the "decay" of the intermediate 
state at rest. 
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Figure 6.1: Left: Experimental data for the photo-nucleon cross-section together with a 
theoretical prediction for the contributions of various subprocesses, from Ref. [17j . 
Right: Comparison of the different contributions to the energy losses of a proton. 



a ~ O.lmbarn and y = 0.5 as 



^GZK ^'^^^^■^^ 400/cm^ X 10^2^ cm^ 2 X lO-^*^ cm~^ (6.3) 



or ^GZK ~ 17 Mpc. Thus the energy loss length of a proton with E > lO'^^ eV is comparable 
to the distance of the closest galaxy clusters and we should see only local sources at these 
energies. Note that it is this smallness of the "horizon scale" that makes the identification of 
UHECR sources feasible, despite the rather poor angular resolution of UHECR experiments, 
5i9 ~ 1°, and deflections of charged particles in magnetic fields. 

The right panel of Fig. 16.11 compares the relative energy losses of a proton due to redshift, 
e~^e~ pair- and pion-production on CMB photons. One clearly recognizes the increase of the 
energy losses by two orders of magnitude around 5 x 10^^ eV and expects therefore a cor- 
responding suppression (not a cutoff) of the cosmic ray flux above this energy. (Remember 
that each shell of thickness Ar contributes the same fraction to the total intensity, if energy 
losses and general relativistic effects can be neglected.) Thus the name GZK cutoff is some- 
what misleading, and the exact strength of the suppression depends on various details as the 
number density of the sources. 

The AGASA experiment detecteded in the 90'ies an excess of events above 10^^ eV com- 
pared to predictions. There had been extensive discussions, if this result requires some kind 
of "new physics", but meanwhile new data of the Hires and Pierre Auger Observatory (PAO) 
confirmed the UHECR flux expected in the presence of the GZK effect. 



Attenuation of nuclei The dominant loss process for nuclei of energy E > 10^^ eV is pho- 
todisintegration ^4 + 7 — > (A — 1) + N in the CMB and the infrared background due to the 
giant dipole resonance. The threshold for this reaction follows from the binding energy per 
nucleon, ~ 10 MeV. Photodisintegration leads to a suppression of the flux of nuclei above an 
energy that varies between 3 x 10^^ eV for He and 8 x 10^^ eV for Fe. 



48 



6.1 Propagation of extragalactic particles 



10 




h180 



Figure 6.2: Deflection map of the GMF for a rigidity of 4x 10^^ V from Ref. [18]; The deflec- 
tion scale at the top is in degree. 

6.1.2 Galactic and extragalactic magnetic fields 

Magnetic fields affect the propagation of charged particles deflecting them, leading in turn 
to synchrotron radiation. While even for ultrahigh energy protons synchrotron losses are 
negligible except in the relatively strong magnetic fields close to sources, the magnitude of 
deflections limits the possibility to perform charged particle astronomy. The deflection angle 
in a regular magnetic field after the distance d is 

. A . 0.52-Z (^) (^) (6.4) 

for a particle with momentum pj_ perpendicular to B. 

Figure 16.21 shows a deflection map for a specific model for the Galactic magnetic field and 
a rigidity of 4 x 10^^ V; for details of the model see Ref. [T8]. The map uses a Hammer- 
Aitoff projection of galactic coordinates with the Galactic center in the middle. The expected 
defiections depend strongly on the direction, but typically trajectories passing the Galactic 
center or plane suffer larger deflections. 

Magnetic flelds beyond the Galactic disk are poorly known and include possibly an extended 
regular and a turbulent fleld in the halo of our Galaxy and a large scale extragalactic magnetic 
field (EGMF). In the latter two cases, the magnetic field is simplest characterized by an r.m.s. 
strength B and a correlation length Ic, i.e. it is assumed that the field is smooth on scales 
below Ic in real space. If we neglect energy loss processes, then the r.m.s. deflection angle 
<5rms = {S'^y^'^ over the distance Ic is 

Rl - W^e\) UoMpcj [iMpcJ [iQ-^ G ) ' ^ ' 

For a calculation of 5i-ms one needs either observational data or reliable theoretical predic- 
tions both for the magnitude and the structure of EGMFs. At present, no single theory for the 
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Figure 6.3: Left: Proton spectrum for different generation spectra [19]. Right: Modification 
factor for the proton spectrum from uniformly distributed sources [20]. 



generation of magnetic field has become widely accepted. Observational evidence for EGMFs 
has been found only in a few galaxy clusters observing their synchrotron radiation halos or 
performing Faraday rotation measurements. The two methods give somewhat different results 
for the field strength in clusters, with B ~ 0.1-1 fiG and B ~ 1-10 ^G, respectively. Outside 
of clusters only upper limits exist for the EGMF. The combination of poor observational data 
and of a missing consistent theoretical picture prevents at present a reliable estimate of the 
infiuence of EGMFs on the propagation of UHECRs. It is likely that future UHECR data 
will provide the first information about EGMF outside the core of galaxy clusters. 
Deflections also imply an average time delay of 



T ~ ~ 1.5 X 10^ 



/ E 



V 1020 eV 



— ) 

lOMpcy 



iMpcy Vio-^G 



B 



(6.6) 



relative to the rectilinear propagation of cosmic rays with the speed of light. As a result, even 
a bursting extragalactic source like a gamma ray burst is seen over time-scales that exceed 
vastly possible observation times. 



6.2 The dip and the Galactic-extragalactic transition 

Transition energy The traditional point of view of the transition between Galactic and 
extragalactic sources assumes that the generation spectrum of both Galactic and extragalactic 
cosmic rays is dN/dEg oc 1/E°' with a = 2.0-2.2 as predicted by Fermi acceleration. In this 
picture, extragalactic sources dominate only above ii^ ~ 3 x 10^^ eV the cosmic ray flux, cf. 
the green line in the left panel of Fig. 16.31 This scenario requires therefore the existence of 
new, isotropically distributed Galactic sources able to accelerate to S ~ 10^^ eV. In a more 
recent view, one uses as generation spectrum of extragalactic cosmic rays dN/dEg oc l/E" 
with a = 2.6-2.7 and can then explain the shape of the cosmic ray spectrum down to a 
fewx lO"*^"^ eV with extragalactic sources [20]. How well the shape is described becomes clearer 
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plotting the "modification factor 77", i.e. the ratio of the diffuse fluxes calculated taking into 
account ah losses and including only redshift losses. In the right panel of Fig. [631 experimental 
data from the HiRes experiment are compared to the theoretically predicted modulation 
factor assuming a simple power-law as generation spectrum. Galactic sources dominate the 
spectrum, when rj > 1: In this case, the observed cosmic ray intensity is larger than expected 
from extragalactic sources alone. Hence Fig. 16.31 motivates the view that Galactic sources 
dominate the spectrum only below a (fewxlO^''-10^'^) eV. 

There exist at least two different explanations for the steep generation spectrum with 
a = 2.6-2.7: Either it indicates that the acceleration to energies above E > 10^'^ eV is not 
mainly based on Fermi shock acceleration, or that the diffuse and the source spectra differ 
not only by propagation effects. In the latter case, it is sufficient that there is a distribution 
of maximal source energies as dn/dEmax oc 1/E^'^. This assumption is plausible, since the 
number of sources able to accelerate to higher and higher maximal energies -Emax should 
decrease with Emax- 

Magnetic Horizon The transition between Galactic and extragalactic dominance in the 
cosmic ray spectrum may be caused by a magnetic horizon hiding extragalactic sources below 
E < 10^^ eV. The maximal distance a cosmic ray can travel is in the diffusion picture given 



where to is the age of the Universe. If we consider cosmic rays with energy below Eq < 10 eV, 
the energy losses are mainly due to the expansion of the Universe, /? = dE/dt = —H. We use 
a "quasi-static" Universe, H{t) = Hq and IoHq = 1. Then E{t) = Eq 6xp(— ffo^) and 



using D{E) = clc/'i{E / EcrY valid for Rl > Ic- (Remember that Ecr was defined by i?L(i?cr) = 



If we assume that a magnetic field with coherence length Ic ~ Mpc and strength S ~ 0.1 nG 
exists in a significant fraction of the Universe, then the size of the magnetic horizon at 
E = 10^^ eV is Thor ~ 100 Mpc. Hence, similar to the GZK suppression above 6 x 10^^ eV, 
we see a smaller and smaller fraction of the Universe for lower and lower energies. As a 
consequence, the spectrum of extragalactic cosmic rays visible to us steepens below 10^^ eV 
and the extragalactic component becomes sub-dominant. 

Energetic reasons suggest that not only the observed diffuse spectrum but also the intrinsic 
spectrum of extragalactic UHECR steepens not far below 10^^ eV: The total luminosity of a 
source depends on its minimal energy as L oc E^^. The combination of large a and small 
-^min cuts across the energy budget of even as powerful and numerous sources as AGNs. This 
problem is softened if there is a transition to an 1 /E"^ generation spectrum at low energies. 

6.3 Extensive air showers initiated by hadrons 

A proton or a nuclei interacting at the top of the atmosphere produces an air shower consisting 
of a core of high energy hadrons that transfer continuously part of their energy to new 



by 




(6.7) 




(6.8) 
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Figure 6.4: Electron number N^. and muon number A^^ (in arbitrary units) for two proton 
and iron showers as function of the depth X. 

electromagnetic sub-cascades, mainly by the process vr*^ — > 27. A smaller fraction of the 
energy is transferred by decays of charged pions to muons and neutrinos. We will not discuss 
these processes in detail, but refer instead only to the results of the simple model discussed 
in Sec. 15. 5t The primary energy Eq, the depth of the shower maximum X^ax and the number 
-^(^max) of particles at the shower maximum are connected by 

N{Xjna.x) oc Eq and X^ax ocln£;o • (6.9) 

Differentiating between primaries According to the superposition model, a nuclei of mass 
number A and with energy Eq interacts as A independent nucleons with energy Eq/A. Thus 
the number of particles at the shower maximum is unchanged but its position X^ax oc IhEq/A 
depends on A. As it is schematically shown in Fig. 16.41 an iron shower develops faster and 
reaches is maximum earlier than a proton shower. Moreover, the shower-to-shower fluctua- 
tions in the case of heavy nuclei as primaries are smaller compared to the ones for primary 
protons. This is very intuitive in the simple superposition picture where a shower initiated 
by a nuclei is described as a superposition of A independent showers. However, the simple 
superposition picture fails to describe quantitatively EAS fluctuations, because of the large 
fluctuations of elementary nucleus-nucleus interactions. 

Figure 16.41 shows also the muon number for a proton and iron shower. First, because 
muons do not multiply through electromagnetic processes (for example, the cross section for 
the process j + A^A + fi'^ + ji^ is strongly suppressed by the muon mass), their number is 
much smaller than the one of electrons and photons. As a result of the slower increase of X^ 
and the large attenuation length, the muon number has a much weaker dependence on the 
depth. This reduces fluctuations and makes, at least in principle, the muon number a better 
estimator for the primary energy. However, this advantage is counterbalanced by the larger 
theoretical uncertainty of and the smaller number of muons. Figure 16.41 illustrates also 
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that observing both and at a fixed depth X provides information about the chemical 
composition of cosmic rays. 

Experimental techniques Present experiments use mainly two different techniques to detect 
ultrahigh energy cosmic rays. The traditional one uses an array of detectors on the ground 
to sample secondaries. Two commonly used detector types are scintillation detectors mainly 
sensitive to electrons and water tanks where the Cerenkov light emitted manly by muons is 
recorded by photomultipliers. The array of detectors is built as sparse as possible, to achieve 
a large aperture within fixed time and money constraints. Therefore only a tiny fraction of all 
particles is measured. Timing information allows one to reconstruct the arrival direction and 
the core of the shower. Then a theoretical formula for the number density p{r) of particles 
as function of the distance to the shower core and of the primary energy is fitted to the 
measurements. Empirically, it is found that p{r) at rather large distance from the core is 
least sensitive to the primary type (i.e. the mass number A of nuclei) and has the smallest 
fluctuation from shower to shower. Depending on the detector type and the spacing of the 
array, this optimal distance varies and amounts e.g. to 1000 m in the case of the PAO. The 
conversion factor between, say p(lOOO), and the primary energy has to be found either from 
Monte Carlo simulations or using different experimental techniques. 

An example for a method that is calorimetric and therefore, at least in principle, allows one 
a model-independent determination of the primary energy is the fluorescence light method. As 
the shower passes through the air, electrons excite nitrogen molecules that emit isotropically 
fluorescence light. Optical cameras on the ground can trace therefore the shower evolution. 
If the excitation probability of nitrogen by a beam of electrons is known, the emitted light 
integrated over the shower trajectory can be used as a measure for the energy of the primary 
particle. Disadvantages of this method are the restriction to moonless nights, the need for a 
careful control of the atmospheric conditions, and the more difficult calculation of the aperture 
of the experiment. 

6.4 Astronomy with charged particles 

Possible anisotropics expected for extragalactic cosmic rays can be classified into four sub- 
classes: 

• At such high energies that deflections in the Galactic and extragalactic magnetic flelds 
are sufficiently small, point sources may reveal themselves as small-scale clusters of 
UHECR arrival directions. This requires additionally a rather low density of UHECR 
sources so that the probability to observe several events of at least a subset of especially 
bright sources is large enough. 

Even if one observes from most sources only one event, a correlation analysis with 
specific sources types may reveal the UHECR sources. This requires however to choose 
a very special class of astronomical sources and/or a high threshold energy, in order to 
keep the number of potential sources within the horizon sufficiently small. 

• Moving to lower energies, the energy-loss horizon of UHECRs and thereby the num- 
ber of sources visible increases. Moreover, deflections in magnetic fields become more 
important. As a result, the identification of single sources is not possible anymore. In- 
stead, anisotropics on medium scales should reflect the inhomogeneous distribution of 
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Figure 6.5: The distribution of normal galaxies color coded according their redshift, in galac- 
tic coordinates. The empty band along the Galactic plane is caused by extinction 
and does not reflect an intrinsic lack of objects |24| . 

UHECR sources that is induced by the observed large-scale structure (LSS) of matter, 
cf. Fig. 16. 5[ Since different source types have different clustering properties, this offers 
an alternative way to identify the sources of UHECRs. 

• At even lower energies, also the LSS of sources disappears, both because the inhomo- 
geneities in the source distribution will be averaged out due to the increased energy-loss 
horizon of UHECRs and because of deflections in the EGMF. Thus the cosmic ray sky 
would appear isotropic, if the Earth would be at rest with respect to the cosmological 
rest frame. As the observation of the CMB dipole shows, this is not case, and a dipole 
anisotropy is expected if the cosmic ray flux is dominated by sources at cosmological 
distance. This is the "cosmological variant" of the Compton-Getting effect discussed 
earlier. 

Since the Sun moves with u = 368 it 2 km/s towards the great attractor, the predicted 
anisotropy is 5cg = (2 -|- 2.7) u ~ 0.6% taking into account the observed spectrum 
I{E) oc E~'^''^ of cosmic rays above the ankle. Deflections of UHECRs in the Galactic 
magnetic field should only displace the dipole axis, but not affect the magnitude of this 
effect. For instance, at energies (2-3) x 10^^ eV and for proton primaries, the dipole 
position should be aligned to the one observed in the CMB within about 10°. Ob- 
serving the Compton-Getting effect at only one energy provides combined information 
on the intervening Galactic magnetic field and the charge of the cosmic ray primaries. 
However, observations at two or more energies break this degeneracy. For example, the 
determination of the average primary charge is straightforward as long as cosmic rays 
propagate in the quasi-ballistic regime and given by the relative shift of the cosmic ray 
and CMB dipole axis at two different energies. 

• Finally, analogous to the latitude effect in the geomagnetic field, the Galactic magnetic 
field (GMF) can induce anisotropies in the observed flux of extragalactic UHECRs, 
for rigidities low enough that blind regions exist. Anisotropies of this kind should be 
expected in models that invoke a dominating extragalactic proton component already 
at £' ~ 4 X 10^'^ eV or extragalactic iron nuclei at E < 10^^ eV [T8]. 

It is not guaranteed that all these four anisotropies can be observed. If EGMFs are large, 
UHECR primaries are nuclei and/or the source density is large, the integrated fiux above 
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Figure 6.6: Chance probability P{S) to observe a larger value of the autocorrelation function 
as function of the angular scale 6 for different combinations of all publicly available 
experimental data; from Ref. [22]. 



the energy where point sources become visible may be for the present generation of UHECR 
experiments too small. Similarly, a transition from galactic to extragalactic sources at a rel- 
atively high energy reduces the chances to observe the cosmological Compton-Getting effect. 



6.4.1 Clustering on medium and small scales 

The cosmological Compton-Getting effect requires that inhomogeneities in the source dis- 
tribution of cosmic rays are averaged out. As the free mean path of cosmic rays decreases 
for increasing energy, anisotropics connected to the large-scale structure (LSS) of cosmic ray 
sources become more prominent and replace the cosmological Compton-Getting effect. The 
exact value of this transition energy depends both on the amount of clustering in the 
source distribution and the free mean path Acr of cosmic rays, i.e. also the primary type. 
For the specific case of proton primaries and a source distribution proportional to the density 
of baryons, Ref. [21] found E'^, ~ 5 x 10^^ eV and a minimal number of order 100 events for 
a detection. 

The first strong evidence for medium-scale clustering was found combining all available 
data sets of UHECR events with energy E > 5 x 10^^ eV and published arrival directions 
until 2005 in [22]. The most convenient tool to search for deviations from anisotropy is in 
this case the cumulative auto-correlation function, 

where N and N' denote the number of pairs with separation angle < -i? in the observed data 
and in a randomly generated data set, respectively. Figure [STUl shows the chance probability 
Pch to observe a larger number of pairs the random data than in the observed data. Clearly, 
the observed arrival directions show a surplus of clustering that is most significant on a typical 
scale of ~ 25°. Similar findings have emerged recently from an analysis of the preliminary 
data from the Pierre Auger Observatory (PAO) [23]. For 64 events with E > 4 x 10^^ eV 
the data presented in Ref. [23] show a surplus of clustering in the broad range from 7 to 
30 degrees. Thus the LSS of matter, probably modified by extragalactic magnetic field, is 
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Figure 6.7: The auto-correlation function w{'d) for galaxies (black), galaxies with M < 24 
(magenta), and AGNs (blue) as function of with la error bars, from Ref. 



reflected in the arrival directions of UHECRs. The obvious question is how one can extract 
the information contained in a plot like Figure \6M 

In order to address this issue, we compare the non-cumulative auto-correlation w{'d), i.e. 
considering the number of pairs N and A^' per bin [d — Af? : + A"!?] , of galaxies and AGNs in 
Fig. 16.71 By construction, w^'d) > indicates overdense regions and one recognizes that the 
clustering at smaller angular scales increases for normal galaxies much slower than for AGNs. 
Tightening the cut in absolute magnitude to M < 24 makes the clustering of normal galaxies 
on the smallest scales comparable in strength to the one of AGNs. Both effects are easily 
understandable: More massive objects are stronger clustered, and the overlap of the brightest 
(optical) galaxies and AGN is large so that these two sets should have quite similar properties. 
A comparison of the auto-correlation function of UHECR arrival directions with the one from 
proposed source catalogues seems to be most promising on intermediate angular scales: First, 
the statistical errors both in cases of sources and cosmic rays become sufficiently small above 
~ 10°. Second, the effect of systematic errors like experimental errors and deflections in 
magnetic fields decreases on larger scales. Turning the second argument around, one can 
use the scale below which deviations between the measured and the expected autocorrelation 
functions appear as measure for the combined effect of deflections in extragalactic and Galactic 
magnetic fields. 

A study of the auto-correlation function of UHECR arrival directions should allow one with 
a rather small event sample to distinguish between GRBs and AGNs as sources. Mainly long 
duration GRBs are discussed as sites for the UHECR acceleration. They occur in star forming 
galaxies that are less clustered as average galaxies. Thus GRBs should be distinguishable from 
AGNs by their clustering properties. 

The strong clustering of AGNs as the paradigm source of UHECRs requires a re- 
interpretation of the traditional view of small-scale clustering. Often, one considered for 
simplicity uniformly distributed sources and distinguished between chance coincidences and 
true clusters, i.e. events from the same source. For clustered sources, the probability to ob- 
serve multiplets from several sources in an overdense region is much increased and this might 
be the main source of small-scale clusters. 
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Figure 6.8: Skymap in galactic coordinates showing the arrival directions of the 27 highest 
energy cosmic rays with E > 5.7 x 10^^ eV (black circles) detected by Auger and 
472 AGN within 75 Mpc (red stars). The blue region indicates the field of view of 
Auger, from Ref. |25j . 

6.4.2 Correlations with astrophysical sources 

The angular resolution of modern cosmic ray experiments reaches ~ 1 degree, and thus even 
for the UHECRs with the highest energies observed, E ~ 10^*^ eV, deflections in magnetic 
fields are likely to be more important. In order to avoid a too large number of potential 
sources per angular search bin, one has thus to choose either a very specific test sample, e.g. 
a small subset of all AGNs, or a very high energy cut such that the cosmic ray horizon scale 
is sufficiently small, say of the order 100 Mpc. Thus we see that the GZK effect is a crucial 
ingredient for a successful correlation analysis. 

Until now, all claims for correlations of UHECR arrival directions with specific source 
classes have remained controversial. It was anticipated that the increased aperture of the 
PAO should improve this situation. The recent analysis of the Auger experiment [25] claims 
indeed correlations between the arrival directions of the highest energy cosmic rays with 
E > 5.7 X 10^^ eV and nearby AGNs. The red stars in Fig. 16.81 mark the positions of AGNs 
with distances < 75 Mpc, while the open circles show the search bin with opening angle 3.1° 
around the CR arrival directions used in this analysis. An overdense region of AGNs along 
the supergalactic plane (dashed line) is clearly visible, where the angular distance between 
AGNs is considerably smaller than the size of the search bin used. The observation of a 
larger number of correlations between UHECR arrival directions and AGNs than expected 
from a random distribution, as in Ref. |25j , should be therefore interpreted only as evidence 
that the UHECR sources have a similar distribution as AGNs. The differentiation among 
different sources (e.g. AGNs vs. GRBs or among different subtypes of AGNs) would require 
a comparison of the correlation (and/or auto-correlation) signal for different source types. 
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7.1 Connecting neutrino, 7-ray and cosmic ray fluxes 

Neutrinos are copiously produced by nuclear reactions in stars and SNe, and as secondaries in 
cosmic ray interactions. Because of the small distance to the production site, solar and atmo- 
spheric neutrinos are experimentally easiest accessible and the only ones regularly detected 
at present. The smallness of their interactions that makes neutrino detection experimen- 
tally so challenging turns neutrinos into ideal messengers to study the interiors of stars, the 
surroundings of supermassive black holes in AGNs or the edge of the universe. 

As we noted in the introduction, any process involving hadronization, i.e. the formation of 
colorless hadrons out of quarks and gluons, leads mainly to the production of pious. Isospin 
symmetry fixes then the ratio of charged and neutral pions. For instance, close to the thresh- 
old, the two possible channels of photo-meson production, p + ^ p + ir^ and p + 7 n + vr^ 
have a ratio that is close to 2:1. While the tt^ decay produces two photons, the final state of the 
second channel contains after the decays n ^ p-\-e~ -\-Ve and vr'*' fi^ + v^ — > + z^e + i^^t + z^e 
two electrons and four neutrinos. Clearly, the production of neutrinos is thus intimately tied 
to the one of photons and electrons, and both depend in turn on the flux of primary cosmic 
rays. 

While photons and electrons start electromagnetic cascades, either on photons in the source 
or from the CMB, neutrinos reach us suffering only from redshift losses. As a result, the 
measurement of X- and 7-rays fluxes can be used to limit possible high energy neutrino 
fluxes. The connection between the observed cosmic rays and neutrinos is more model- 
dependent. A popular set of assumptions is the following one: Protons are accelerated in a 
source region, confined by magnetic fields. They interact on background photons (say in the 
UV range) in the source, neutrons from p + 7uv — > n + vr"*" escape and provide the observed 
cosmic rays, while the Tr"*" decays produce neutrinos. This argument does not exclude that 
(a subset of) sources exists with sufficiently large interaction depth r such that also neutrons 
interact frequently inside the source. The cosmic ray flux from these sources would be strongly 
suppressed and therefore this type of sources has been coined "hidden sources." Before we 
discuss the connection of neutrinos, cosmic rays, and 7-ray fluxes, we consider the yield of 
neutrinos and photons for a single source. 

Photon and neutrino yields at the source The yield Yi{E) of secondaries is defined as the 
dimensionless ratio of the secondary flux (f)i{E), i = 7, z^, and the product of the initial proton 
flux (j)p{E) and the interaction depth r, 



We consider only in detail the simpler case of photons as secondaries. Generalizing Eq. (jS.lOp . 
the total photon flux from neutral pion decays that are in turn produced as secondaries in pp 



y^{E) 



r(t>p{E) ■ 
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or p7 scatterings is 



dE,NAEr,)'^^%^r{E,,^, 



pi 1 



(7.2) 



where A'^o denotes the multiphcity of neutral pions and r the interaction depth. For the 
same assumptions as used in the derivation of Eq. (j5.17|) . we obtain as approximation for the 
photon yields 



Y,{E) 



a 



(7.3) 



where now in the definition of Zp^o (a) the neutral pion multiplicity is included. The calcula- 
tion of the neutrino yields proceeds along the same lines. However, the resulting expressions 
are rather lengthy, because the decays are non-isotropic, the decay chains are more complex 
and the final-state particles massive. The relevant expressions can be found in Refs. [E1,E2]. 

In Fig. 17.11 the neutrino yields Yy{E) from various subprocesses (pion, kaon and charm 
decays) are shown for a thermal distribution of photons with temperature T = 10^ K and 
three different interaction depths. As long as r < 1, i.e. multiple scatterings can be neglected, 
the yields are by construction independent from r. The yields are also independent from T, if 
one uses instead of the energy E the parameter x = Eto/mp, where uj = l.QT is the maximum 
of the Planck distribution. Therefore they are a convenient tool to estimate photon or neutrino 
yields as long as r < 1. Tabulated values of Y for different a can be found in Ref. [E1,E2]. 




Figure 7.1: Neutrino yield Y as function of x 

and neutral kaon, and charm decays for r 



Euj/mi for T 



.p ^ — 10 K from pion, charged 
0.05,0.5 and 50; from Ref. 



Cascade limit The general concept of electromagnetic cascades discussed in Sec. 15.51 has 
several important applications. Firstly, we can use it to derive a simple bound on the UHE 
neutrino fiux. Using the branching ratio 2:1 for p + ^ ^ p + tt^ and p + ^ ^ n + 7r+, the 
approximation that the four leptons produced in muon decay share equally the energy as well 
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as the fact that in neutron decay almost ah energy is taken by the proton, the ratio of energy 
transferred to neutrinos and the electromagnetic channel follows as 1:3, 

If we approximate the neutrino intensity in the range [i?, E'max] by a power-law, luiE) = 
AE~"-, then the integral neutrino intensity is connected to ujy via 



Iu{> E) " 



a + 1 



1 



E 



-1 



(74) 

47r E ' ^ ^ ^ 



where LOy is the energy density in neutrinos in the same energy range \E^E^^-^. The first 
two terms are bounded by one, provided that the energy range considered is large enough, 
Q;ln(£'/£^max) > 1- With ojy < Sco^y, we obtain thus as bound for the integrated neutrino flux 

«>E)<^^. (7.5) 

In order to apply this bound, we have to distinguish between different possible photon 
backgrounds on which cascades can develop. The only universal one is the CMB leading 
to a guaranteed flu^lll of UHE neutrinos via the GZK reaction of UHECRs as well as to a 
contribution to the diffuse photon flux in the 100 MeV-10 GeV range. This energy range 
has been observed by the EGRET experiment and thus a bound on the diffuse extragalactic 
gamma-ray flux can be derived. Integrating this maximal flux, the limit on uj^ follows as 

cj^ < 2 X lO^'^eV/cm^ (7.6) 

If we assume that Iy{Ey) = AE~'^, we can combine Eqs. ()7.5p and ()7.6p into a limit on the 
differential neutrino intensity. 

El Iy{Ey) < 2 X 10"*^ GeV cm'^ sr"^ . (7.7) 

Electromagnetic cascades may happen also inside many sources. Since the Thomson cross 
section is a factor 3000 larger than the cross section for photo-meson production, (TTh/cp7 ~ 
3000, sources where photo-meson reactions are relevant are likely optically thick. To be 
specific, let us assume that the source of opacity are UV photons, ~ 10 eV, as it is the case 
for many AGNs. Hence the critical energy for pair production is reduced to E^r ~ 25 GeV 
compared to E'er ~ 400 TeV for cascades on CMB photons. Thus the picture of the escaping 
photon radiation is as follows: 

• Photons with E ^ E^r escape from the source only from an outer layer of such a 
thickness I that the interaction depth r^^ = la^^yU^ for pair-production is small, r^^ ~ 1. 

• The source is bright in the ~ 10 — 100 GeV range, where the last generation of photons 
generated in the cascade can escape. 

• Electrons with energy below -E'er continue to scatter in the source in the Thomson 
regime, producing photons in the X-ray range with average energy 

P = ~ lOOeV-O.lMeV (7.8) 

3 

for Eg = 1 GeV and 25 GeV, respectively. 



^If UHECRs are dominantly nuclei, a neutrino fiux at lower energies from neutron decays (via A + 'ysK 
{A~ l)+n) results. 
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Cosmic ray limits The total energy density in neutrinos can be directly bounded by the 
energy density injected in cosmic rays by extragalactic sources, if neutrinos are produced 
only in transparent sources. If one assumes further that only neutrons can leave the sources, 
this bound becomes a prediction. To proceed, one has to specify the generation spectrum of 
extragalactic cosmic ray, i.e. one has to choose essentially between the traditional point of 
view (injection spectrum dE/dN oc E~'^ and transition energy around 3 x 10^^ eV) or the 
dip scenario [(effective) injection spectrum dE/dN oc E~'^'^ and transition energy around 
10^^ eV). Using the first choice, several authors advocated bounds or expectations for the 
neutrino flux of order 

El I^{E^) ~ few X lO^^GeV cm^^ g-i g^-i _ ^7 9) 

Thus the cosmic ray limit is about a factor 100 stronger than the cascade limit. However, one 
should be aware that even within the framework of transparent sources cosmic ray bounds 
obtained for a fixed set of parameters can only be indicative, because the sensitivity of the 
cosmic ray and the neutrino fluxes on several free parameters like the redshift evolution of 
their sources differs strongly. 




Figure 7.2: A collection of various bounds on neutrino fluxes: The solid line at the top cor- 
responds to the cascade limit (17. 7p . the lower dashed horizontal line to the cosmic 
ray limit (|7.9p . while the lines in-between correspond to varying assumptions used 
for the cosmic ray limit, from Ref. [27j . 



7.2 Neutrino interactions, masses and mixing 

One of the most surprising properties of the weak interactions is that only left-chiral fermions 
{P-ip) and right-chiral anti-fermions {ipPj^) participate in charged-current (CC) weak inter- 
actions, 

CcC = E h,aryL,aW- + h.C. . (7.10) 
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Thus weak interaction violate parity P leading to a doublet /singlet structure of fermions, 

, - , , , ^ , eR,^lR,TR. (7.11) 

/L \^ /L 

The current ()7.10p defines what we mean with electron neutrino: It is the neutrino created 
together with an electron in a CC interaction. 
Since ordinary Dirac mass terms 

m/VvV'/ = "i/(V'/V/ +V'/V'/) (7.12) 

are not gauge invariant, one introduces fermion masses in the standard model by gauge- 
invariant Yukawa interactions. These interactions need not to know about weak interactions 
and are thus not diagonal in flavor space, 

jOrrn^s = - E ^L,aMi"^Un^f, - ^ I L ,o.M^^p'^ I R^p + h.C, (7.13) 

i.e. the mass matrices Mq,^ are not hermitian or even diagonal in the basis of the weak eigen- 
states. (We denote the weak eigenstates Va = {i^e, i^r} and la = {e, r} by Greek indices 
a, /3 . . ., and the mass eigenstates by Latin indices . . . Furthermore, we have assumed that 
the neutrinos have only Dirac mass terms to simplify the formulas). Since the mass matrices 
are not hermitian, they cannot be diagonalized by a simple unitary transformation. However, 
arbitrary mass matrices can be diagonalized by a biunitary transformation, 

Mj^-l^ = [7('^)tM(^)r('^) , (7.14) 

where [7('^)f/('^)t = 7"{'^)7"{'^)t = \ Then, the connection between weak and mass eigenstates 
is given by 

l'R,a = J2 ' ''L,a = J2 ^at ^L,i , (7.15) 

i i 

and similar equations hold for the charged leptons. 

Inserting the transformations of Eq. (j7.15p into the charged current Lagrangian (j7.10p of 
the SM, results in 

Ccc = -4nH h,a''Vi,UL,jW- + h.C, (7.16) 

where V = [/('^^)t[/('^) is the analogue of the CKM matrix in the lepton sector. Since the 
charged current interaction involves only left-chiral fields of both charged leptons and neu- 
trinos, the product of the two mixing matrices of the right-handed leptons, T^^^^'^T^'^\ is 
unobservable. The mixing V = JJ^^^^'^U^'^^ between charged leptons and neutrinos has two 
different sources: It could be ascribed either completely to mixing in the neutrino {V = U^'^'^) 
or in the charged lepton sector {V = U^^^'^^), or most probably to some superposition of both. 
Only the combination V is observable and, by convention, we set urn = 1. 

In the case of massless neutrinos we can choose the neutrino mass eigenstates arbitrarily. 
In particular, we can set any given U^^^\ hereby rotating away the mixing. 

This shows that neutrino masses are a necessary condition for non-trivial consequences of 
mixing in the lepton sector. 
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Figure 7.3: Production of a superposition of neutrino mass eigenstates vi in pion decay and 
subsequent detection of the neutrino flavour via the secondary lepton l'^. 



Neutrino oscillations in vacuum 



Let us consider e.g. neutrinos produced in charged pion decay. The ratio of tt ^ efg and 
TT — > ^^v^ decay rates is 



R 



r(7r ei^e) mg (m^ — mr^ 



2\2 



r(7r ^ IIU^) 777,2 (^2 _ ^2 



2^2 



1.28 X 10 



-4 



(7.17) 



since angular momentum conservation in the pion rest frame requires a helicity flip of the 
lepton. Similar, in neutron decay and in fusion reactions in stars only are emitted, because 
of energetic reasons. Hence, in many occasions we start with a (nearly) pure flavor state. 

The time-evolution between creation of an arbitrary state at t = and detection at t 
becomes simplest, if we decompose the weak interaction eigenstate into mass eigenstates 

IKt)) = E^mV.)e-'''*- (7.18) 



Neutrinos are in all applications ultra-relativistic, 

Ei = (/ + rnf )i/2 ^ p ^ m|/(2jo) , 

where we have assigned also a definite momentum to the states 



(7.19) 



(7.20) 



The probability for a transition from the flavor to after the distance L = rf is 



J2 Upkew{-iEt)Uak 

k=l 



(7.21) 



where we introduced also Am?- = |m? — m3\ 
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To be specific, we consider the simple case of two-flavor neutrino oscillations. Then the 
survival probability is 

P^,^^, = 1 - sin^ fir sin^ 2^ (7.22) 



with the vacuum oscillation length 



E f E \ I e\ 

lose = 47r- — 5" K, 2.5 m 



Thus the survival probability depends only on L/E, the mass squared difference of the two 
neutrinos, and one mixing matrix element. Most strikingly, the smallness of the neutrino 
masses makes it possible to observe the interference of quantum states on macroscopic length 
scales. 



7.3 Atmospheric neutrino oscillations 



The principle production modes for neutrinos produced by comic rays in the atmosphere are 

y^(t>^) and ^i^ ^ + i'e{i>e) + v^{v^), 



pion, kaon and muon decays. With tt^ ^ /i^ + i^^{i^fj.) and fi^ —>■ + t'e(fe) + i^fj(i^^t), one 



expects the ratio 

I^e + i^e Ne 1 



(7.24) 



if all particles can decay. Muon decays become suppressed for E > 2.5 GeV (A^ > Xatm) 
and thus the ratio should decrease for higher energies. Since most neutrinos are produced 
by low energy cosmic rays, geomagnetic effect must be included in the flux calculations. On 
the experimental side, the connection of the visible energy and direction (of the secondary 
electron or muon) to the one of the neutrino can be made only on a statistical basis and 
requires good knowledge of detector response and acceptance. 

Part of these problems can be avoided by looking at the double ratio 

data /^7 nr\ 

^ " (iVc/A^,)MC • ^^^ ^^ 

If in the Monte Carlo no oscillations are included, one expect i? ~ 1 if neutrinos are massless. 
The first conclusive deviation from i? = 1 was found by the Superkamiokande experiment in 
1998 with R = 0.69 ± 0.06. 



Ratio of up- and down-going muons A decisive test for the oscillation hypothesis is the 
ratio of up- and down-going muons: The oscillation length varies between 20km and 12800km, 
and as function of the zenith angle the ratio should show the form predicted by the neutrino 
oscillation formula. By fitting the experimental data to the free parameters Am^ and sin^ i? 
of the neutrino oscillation formula, one derived Am? ^ 2.2 x 10~^ eV^ and sin^ ~ 1. 

The atmospheric anomaly was later confirmed by sending a beam from KEK to Kamioka. 
The distance KEK-Kamioka is L = 250 km, while the mean neutrino energy was P ~ 1.3 GeV. 
The combined result of the two experiments is shown in Fig. 17. 4^ yielding the surprising result 
that the atmospheric neutrino mixing angle is - in contrast to the quark sector - close to 
maximal. 
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Figure 7.4: Overlay of the Superkamiokande and the K2K result. 



7.4 High energy neutrino scattering 

Parton model and scaling violation Neutrino-nucleon scattering can be described in the 
parton model. This model assumes that a nucleon probed at sufficiently high > Qmin ~ 
1 GeV^ behaves as a collection of independently interacting partons, i.e. quarks and gluons. If 
the function fi{x) denotes the probability to find a parton of type i with momentum fraction 
Pi = xpN and a{s) the corresponding parton cross section, then the neutrino-nucleon cross 
section follows as ^ 

a{s)=J2f dxMx)a{s). (7.26) 



In its original formulation, the parton model assumed that the parton distribution functions 
fi{x) depend only on x. We can convince ourselves that fi{x) should depend also on by 
looking at a Bremsstrahlung process. Consider a quark as an external line of a Feynman 
diagram with momentum p and arbitrary mass m emitting a gluon of momentum q, 



{p + qY — fn^ '^pq 2ujE{1 — cos t?) ' 

where oj is the energy of the gluon, E and v the velocity of the parton and "d the emission 
angle. There is a soft divergence for all m, if the energy of the emitted gluon goes to zero, 
tj ^ 0. Additionally, there is a collinear singularity for light partons, f ^ 1, when — > 0. 
These singularities generate logarithmic divergences integrating over the final phase space, 

d(7„+i = ln2(QVQmin)da„ , (7.28) 

that can compensate the strong coupling although is small for large . As a result, parton 
cascades or jets develop in fragmentation processes. Similarly, the parton number inside a 
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Ej, [GeV] 

Figure 7.5: Cross sections for uN interactions at high energies: dotted hne, a{vN — > + 
anything); thin hne, (t[i'N — > ^~ + anything); thick hne, total (charged-current 
plus neutral-current) cross section, from Ref. |28j . 



nucleon increases with increasing . Although the neutrino-quark cross section is constant 
above s » rn^f, the neutrino- nucleon cross section continues to grow as c oc E^-"^, because 



the number of accessible partons with x > 



Qmin/^ increases as xq{x,Q ) ~ x 



-0.33 



DGLAP equations The Q^-dependence of the parton distribution functions is described in 
perturbative QCD by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations 



dqi{x,Q^) 
dt 

dgiix,Q'^) 
dt 



asm 

27r 



y 
y 



+ g{y,Q'')P, 



/ X 

''\y 

X 

y 



(7.29) 
(7.30) 



with t = ln((5^/(5o) Is = J2]jLi[Qj + Qj]- The splitting functions Pij{^), with i,j = q,g 
give the probability that parton j with momentum y radiates a quark or gluon and becomes 
a parton of type i with fraction (|) of the momentum of parton j. 

If we compare this equation to the production term of the transport equation for cosmic 
rays, Eq. (j3.38p . we see that the structure of the two equations is the same. The reason for 
this formal agreement (or the simplicity of the DGLAP equations) is that interference terms 
can be neglected in the QCD cascade, since it corresponds to a semi-classical evolution. 



7.5 Astrophysical sources of high energy neutrinos 

Experiments and techniques Main problem of high energy neutrinos physics is the atmo- 
spheric neutrino background and the expected low fluxes from astrophysical neutrino sources. 
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Figure 7.6: Averaged shadow factor S{E) as function of the neutrino energy. 



The search for neutrino sources is most promising at high energies, because the atmospheric 
neutrino flux is very steep: At high energies, pions scatter and muons do not decay. Main 
detection methods are at present Cherenkov light (in ice or water) emitted e.g. by a relativis- 
tic muon produced in the reaction + — > /u + X. In the future, the radio or acoustic 
signal produced by a horizontal shower initiated by a neutrino reaction may allow one the 
supervision of even larger detector volumes. 

The event rate of muons with energy E > ii^min in a detector of area A is 



/ 



Rate = A / dEdO P^{E; E'^''')S{E) I^{E^ 



(7.31) 



Here, the shadow factor S takes the attenuation in the Earth into account, and Pf^iE; E^ 
is the probability that a muon is created and reaches the detector with > £'min- 



ti ) 



Attenuation and Shadow factor The growth of the neutrino-nucleon cross section with 
energy means that above Ei, ~ 40TeV the Earth becomes opaque even for neutrinos. If 
we neglect any regeneration effects and assume that the flux is isotropic, this attenuation 
can be represented by a shadow factor. The shadow factor averaged over the upward-going 
hemisphere is 



1 /"O 

and is shown in Fig. 17.61 



dcosT? / d<p exp[-A(i9)/Aint(E)] 



(7.32) 



Average range of muons The second ingredient in the rate calculation is the probability 
that the neutrino creates a muon that is energetic enough to arrive at the detector with an 
energy E^ larger than the detector's threshold energy ii^™™. The probability that a muon can 
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[GeV] 



Figure 7.7: Average ranges (in rock) for muons produced in charged-current interactions of 
neutrinos with energy Ei,, at threshold energies -E"™™ = 1 and lOTeV. 



be recorded in a detector depends on the average range (R) of a muon in rock 

{R{E,- Ef^)) = \— I dy R{E,, E^n^ ^^CciEu.y) ^ ^^ 33^ 

^ acc{Eu) J ^ dy 

where the energy of muons produced in a charged-current interaction of neutrinos with matter 
is E^ = (1 — y)Ey. After a high energy muon is produced, it undergoes a continuous energy 
loss as it propagates. The range R of an energetic muon follows from the energy-loss relation 

- AEjdX = a{E^) + /3{E^)E^, (7.34) 

where X is the thickness of matter traversed by the muon in g/cm^. The first term represents 
ionization losses, while the second term represents bremsstrahlung, e~^e~ pair production and 
nuclear interactions. If the coefficients a and /3 are independent of energy, we can approximate 
their values to be a = 2.0 x 10~^ GeVcmwe"^ (cmwe = g/cm^) and /3 = 3.9 x 10~® cmwe~^. 
Integrating Eq. (j7.34p . the muon range is 

R{E., Sf") = ^(^f - X{E.) = yln . (7.35) 

The average range of muons from charged current neutrino interactions is shown in Fig. 17.71 
for threshold energies 1 TeV and 10 TeV. A muon produced with E^ = 10 TeV will travel on 
average a few kilometers until its energy is degraded to 1 TeV. 

The probability that a muon neutrino of energy E produces an observable muon is then 

P,iE, Ef^) = Na a{E){RiE; Ef^) (7.36) 

and is shown in Fig. 17. 8i A useful approximation for energies E > TeV is P^(E,E^^'^) ~ 
2 X 10-^{E^/TeYf-^. 
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Figure 7.8: Probability that a neutrino of energy produces an observable muon with energy 
exceeding = ITeV and lOTeV. 



Sensitivity and event rates We derive now a rough estimate of the event rate expected 
in a neutrino telescope. We assume -Emin = 10^ GeV as minimal energy in order to avoid 
the background of atmospheric neutrinos and a neutrino intensity equal to the cascade limit 
(17. 7p . Then we approximate the remaining ingredients of the integrand in Eq. (I7.3ip as 
~ 10-3(^/£;^i„)0-8 and S ~ {E / E^^i^)-^-^ . Choosing as area A = (0.1km)2, i.e. the 
area covered by the current ICECUBE prototype AMANDA, we obtain of order 1 event/year. 
Thus the sensitivity of this experiment corresponds roughly to the maximal flux allowed by 
the cascade limit, while an 1 km'^ neutrino telescope like ICECUBE will probe fluxes down to 
the ones compatible with the cosmic ray limit. 



7.6 Exercises 



1. Derive from Fig. 6.3 the energy density of (extragalactic) cosmic rays above 10^^ eV for 
a = 2 and 2.7. Estimate from Fig. 6.1 the "life-time" of a cosmic ray and thus the 
required emissivity C. 

2. Derive for two neutrino the oscillation probability, Py^^y^, starting from the mixing 
matrix 

TT - ( sint?e"^ \ 

~[ -sint? cos^e^f' J ' ^'■'^'> 

Show that that the Majorana phase (j) does not enter the oscillation probability. Con- 
vince yourself that probability is conserved. 
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8.1 Air shower and new particle physics 

At present accelerator data from the Tevatron constrain new physics directly only up 1 TeV, 
while in the near future LHC extends this range up to 14 TeV. Therefore one might wonder 
if cosmic ray data can be used to constrain new particles physics at energy scales that are 
not accessible to human accelerators. However we know from the discussion in Sec. 2.2 that a 
significant change of hadronic cross sections requires the modification of scattering with small 
four-momentum transfer t - while new physics is most often connected to the regime of large 
masses and therefore large t. More promising is neutrino scattering where at high energies 
the whole range < \t\ < m^r contributes uniformly to the neutrino-nucleon cross section on 
the parton level. 

IMew resonances In sees. 2.2 and 7.5 we considered mainly the i-channel contribution to 
cross sections. To complement this discussion, we assume now that a new particle with mass 
m and decay width T acts as a s-channel resonance in fN scattering. (Concrete examples 
are models containing lepto-quarks or a sneutrino in i?-parity violating super symmetry.) In 
the narrow-width approximation, F <C m, for the Breit-Wigner denominator 

1 2 

(s — m^)^ + m^r^ mT ^ ^ 

in the partonic cross section, the integration over x just picks out partons with s = xs = rri^. 
The total cross section is 

(Ti,n{s) ^^—;^xqi{x = /s,Q'^ = s) , (8.2) 

and hence is suppressed for large m. Experimental limits require either a small coupling g, 
a large mass m, w? > or a suitable combination. This simple example illustrates that 
it is difficult to increase the interactions of neutrinos by a large factor. An exception may 
be models with large extra dimensions. An important ingredient for the large increase of the 
uN cross section in these models is that the neutrino couples not only to quarks but also to 
gluons. 

Neutrinos as the particles with the smallest cross section are the most sensitive tool to 
search with cosmic ray data for deviations from the standard model predictions caused by 

unknown physics. They produce mostly showers close above or below the horizon. Looking 
at the zenith angle dependence of the event rate, one can disentangle the cross section and 
flux, and thus test for deviations from the standard model prediction. 

Understanding strong interactions We have seen that several aspects of strong interaction 
cannot be calculated within QCD but are described by phenomenological models. Moreover, 
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Figure 8.1: Experimental data (red, with errorbars) and model predictions for Xmax distri- 
bution as function of energy, from Ref . [29j . 

even usual methods and results from perturbative QCD may become unreliable at ultrahigh 
energies. A simple example are the DGLAP equation and the parton distribution functions: 
In the standard picture, only splittings X — > YZ are taken into account, and for decreasing x 
the number of partons in nucleons increases monotonically. As a result, the neutrino nucleon 
cross section increase like a power, while for s — > oo only a logarithmic growth is allowed. 
Physically, the density of partons cannot become infinite and at some point recombinations, 
YZ X, will balance splittings. 

Although the physical picture seems to be rather simple, its formalization is highly non- 
trivial and has led e.g. to ideas like the "color glass condensate." A test of one of these 
proposals against experimental cosmic ray data is shown in Fig. 18.11 The two red curves 
are the prediction of an often used model for hadronic interaction in cosmic ray physics for 
protons (upper) and iron (lower line). These two curves enclose the experimental data (red, 
with errorbars) and thus this model predicts that the cosmic ray flux is a mixture of protons 
and heavier nuclei. The other two lines are models where the growth of the parton density is 
slower than in perturbative QCD. Both lines are for protons and the corresponding lines for 
iron would be even lower. Hence the model producing the magenta line is already excluded 
by cosmic ray data, while the other one predicts a pure proton flux. 

In summary, cosmic ray data can be used at present only to exclude extreme model predic- 
tions, since the chemical composition of cosmic rays is not independently fixed. The discovery 
of cosmic ray sources and the subsequent study of cosmic ray deflections in magnetic fields of- 
fers an alternative way to measure the chemical composition, providing in turn the possibility 
to test QCD models more precisely at ultrahigh energies. 

8.2 Relics from the early universe as cosmic rays sources 

Evidence for dark matter comes from flat rotation curves of galaxies, the virial mass of 
galaxy clusters, big-bang nucleosynthesis, structure formation and the CMB: Combining all 
these pieces of evidence, one concludes that there exists around six times as much cold dark 
matter (CDM) than baryons. Many well-motivated dark matter candidates that can have the 
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Figure 8.2: A selection of dark matter candidates in the plane cross section versus mass; blue 
and black corresponds to thermal and non-thermal as main production channel. 

correct abundance have been proposed, with a rather big variety in their properties as shown 
in Fig. 18.21 Therefore most of them have to be searched for by some specialized methods. We 
shall discuss only two candidates in more detail that can serve also a source of cosmic rays. One 
of them are "WIMPs," i.e. particles with masses and interactions similar to the weak scale. 
The second one are superheavy DM particles (SHDM) with masses around 10^^ GeV. Before, 
we comment briefly om the other candidates shown in Fig. 18. 2t Neutrinos with interactions 
as predicted by the standard model would have either with J2i=i ~ 10 or few GeV 
the correct abundances as thermal relics. Both possibilities are excluded by various reasons. 
Axions were proposed as solution to the strong CP problem. They are pseudo Goldstone 
bosons that mix with pions, and therefore their masses and decay constants are connected 
by ma = mTjfT^/fa- Gravitinos and axinos are supersymmetric partners of the axion and the 
graviton, respectively. They can be the lightest supersymmetric particle that is protected in 
many models by a symmetry and thus stable. Last but not least, the lightest neutralino, a 
mixture of the supersymmetric partners of higgses, Z and photon, is a typical example for a 
WIMP. 



8.2.1 WIMPs as thermal relics 

When the number A'^ = nV of a particle species is not changed by interactions, then the 
expansion of the Universe dilutes their number density as n oc i?^^. The corresponding 
change in time is connected with the expansion rate of the universe, the Hubble parameter 
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Figure 8.3: Illustration of the freeze-out process. The quantity Y = nx/s is nx divided by 
the entropy density s oc to scale out the trivial effect of expansion. 



-3Hn. 



(8.3) 



H = R/R, as 

dn dn dR ^ R 
dt ~ dR~dt ~~ ^R 

Additionally, there might be production and annihilation processes. While the annihilation 
rate Pv? = (aannv) has to be proportional to v?, we allow for an arbitrary function as 
production rate ip, 

-3Hn - (ir^ + il). (8.4) 



dn 
~dt 



In a static Universe, dn/dt = defines equilibrium distributions rigq. Moreover, detailed 
balance requires /3ng„ = and thus we can eliminate the unknown function ip, 



dn 
dt 



— = -3Hn - (dann?^) ("-^ 



n, 



eq 



(8.5) 



This equation together with the initial condition n ~ rieq for T ^ oo determines n{t) for a 
given annihilation cross section cJann- The evolution of nx (divided by the entropy density 
s oc 1 /R^ to scale out the trivial effect of expansion) is shown schematically in Fig. 18.31 As the 
universe expands and cools down, nx decreases at least as Therefore, the annihilation 

rate oc quenches and the abundance "freezes-out:" The reaction rates are not longer 
sufficient to keep the particle in equilibrium and the ratio nx/s stays constant. 
Numerically, one obtains for the relative abundance ^Ix = Px/Pcr of CDM 



Qxh^ 



2 X 10-^^cm7s 



■ Xf , 



(8.6) 



73 



8 Cosmic rays as tool for particle physics 



where Xf is the ratio of the mass and the freeze-out temperature, Xf = M/T{. Thus the 
abundance fix is inverse proportionally to the thermally averaged annihilation cross section 
(ft'ann), since a more strongly interacting particle stays longer in equilibrium. The abundance 
depends only logarithmically on the mass m via the freeze-out temperature Xf . Typical values 
of Xf found numerically are Xf ~ 20. Prom the observed value Ocdm ~ 0.20 one obtains 
the annihilation cross section in the early universe — that controls also the annihilation rate 
important for indirect searches today. The size of the cross section corresponds roughly to 
weak interactions and has led to the paradigm of "WIMPs" , weakly interacting massive dark 
matter. 

Unitarity bounds each partial-wave I of the thermally averaged annihilation cross section 
as {(Tajinv)^''^ < const. / {vm^) . Requiring $1 < 0.3 leads to m < (20-50) TeV. This bounds the 
mass of any stable particle that was once in thermal equilibrium. 

8.2.2 SHDM as inflationary relic 

The energy density of non-relativistic particles decreases only as px oc while the one of 

radiation decreases as p oc Therefore particles that never came in chemical equilibrium 

with radiation should be created only in tiny amounts in the early universe. Already the 
production by gravitational interactions at the end of inflation is suflicient to ensure an 
abundance Clx ~ 1 for stable particles with masses mx ~ lO^^GeV. 

As an illustration, we consider the Klein- Gordon equation for the field modes 0^ of a scalar 
field in the Priedman-Robertson- Walker metric, 

Mr]) + mlsirj)Mri) = , (8.7) 
where the eff'ective mass depends on an unknown parameter ^ and is 

mlsiv) = k^ + MlR^ + (6^ - 1)| . (8.8) 

Since rriefr is time dependent, vacuum fluctuations will be transformed into real particles. 

Thus, the expansion of the Universe leads to particle production. 

The predicted abundance of X particles today from inflationary cosmology is 

where Trm is the temperature at reheating. Since the production is gravitational, the result 
is modcl-dcpcndcnt only through the dependence on cosmological parameters like the value 
of Hubble parameter at the end of inflation and the reheating temperature Trh- 

8.3 Indirect detection of dark matter 
8.3.1 Dark matter in the Galactic halo 

The average density of DM in the Galaxy is strongly increased compared to the extragalactic 
space, nyiw/n-cx ^ 10^. Therefore the annihilation rate of DM can become again appreciable 
inside the Milky Way, and in particular, in objects where DM is strongly accumulated. 
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DM distribution and clumps Typical results from AT-body simulations are for the smooth 
DM mass density profile psm 

with pq = 0.3GeV/cm^ as the dark matter density at the solar distance from the Galactic 
center, and a ~ 25kpc as the characteristic scale where the slope changes from r~"^ to 
( "Navarro- Frenk- White profile"). At small radii r < 1 kpc, the missing resolution of A'^- 
body simulations, the influence of baryonic matter and of the galactic SMBH make a reliable 
estimate of the DM density difficult. 

According to the model of hierarchical structure formation, the first objects to form were 
the smallest structures. For the case of WIMPs, 10^^ Earth-mass DM halos about as large 
as the Solar System might be in the Milky Way. An important contribution to the total 
annihilation signal can be given by small clumps that have a denser core and may be at small 
distance. 

Photons, neutrinos and antiprotons from DM annihilation in the halo The secondaries of 

DM annihilations will be the stable particle of the standard model, i.e. photons, neutrinos, 
electrons and protons. For the latter two, only the anti-particles may provide some useful 
information. 

The differential flux of the final state i at the Earth from DM annihilations is 

iUE,^) = ^pif d,^^Mll, (8.11) 

dE 2 nix Ji-o.s. 47r 

where r(s, tjj) = (r,?) + — 2 Tq s cos ipY^"^ , ip is the angle between the direction in the sky and 
the galactic center (GC), Vq « 8.0 kpc is the solar distance from the GC, and s the distance 
from the Sun along the line-of-sight (l.o.s.). In terms of galactic latitude b and longitude I, 
one has cosV' = cos 6 cos L The energy spectrum dN^/dE can be calculated only within a 
specific model for the DM particle. 

SHDM as UHECR source The original motivation to introduce SHDM was the "AGASA 
excess" , i.e. a surprisingly large number of cosmic ray events with energy > lO^'^ eV and 
thus above the GZK cutoff. Since SHDM behaves by definition as CDM, its abundance in 
the galactic halo is strongly enhanced. As a result the cosmic flux is dominated by the halo 
component and the GZK cutoff is absent in this model. 

The hadronization spectra of superheavy particles can be reliably calculated using standard 
QCD methods. The predicted spectrum in the SHDM model, dN/dE cx E~^-^, cannot fit 
the observed UHECR spectrum at energies E < (6-8) x 10^^ eV. Thus mainly events at 
E > (6-8) X 10^^ eV, and most notably any excess at energies beyond the GZK cutoff, could 
be produced by SHDM decays. As for all hadronization processes, the main component of the 
UHE flux are neutrinos and photons from pion decay. As additional signature, one expects 
a Galactic anisotropy, because the Sun is not in the center of the Galaxy. The degree of 
this anisotropy depends on how strong the CDM is concentrated near the Galactic center 
- a question under debate. First results from the PAO find neither the predicted Galactic 
anisotropy nor the dominance of photon primaries at the highest energies. If confirmed, 
SHDM can play only a sub-dominant role as source of UHECRs. 
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8.3.2 Neutrinos from dark matter annihilations in the Sun and the Earth 

Dark matter particles X scatter on matter in the Sun or Earth, lose energy, and may become 
gravitationally bound. They continue to scatter, gaining or losing energy. If energy looses 
dominate, they sink down to the center where they eventually annihilate. The directed signal 
of high energy neutrino would provide a rather clear evidence of DM. In the following, we 
illustrate the basic steps in the calculation of the resulting neutrino flux, following rather 
closely Ref. [30]. 

The total number Nx of dark matter particles in a celestial body is determined by three 
processes: The capture rate Rc, the evaporation rate Re, and the annihilation rate Ra, 

Nx = Rc- ReNx - RaNI = Rc-Te-Ta. (8.12) 

If the DM particles are self-conjugated, X = X, each annihilation reduces their number by two 
and thus one should replace RaN\ by 2RaN\. The evaporation rate takes into account that 
an already captured dark matter particle might be "kicked out" in a scattering process and is 
only important for particles not much heavier than nucleons, say mx ^ fewGeV. Therefore 
we will set Re = and obtain as solution of Eq. (j8.12p 

Nx{t) = iVeqtanh(t/r) (8.13) 

with r = {RcRa)~^^'^ as time-scale to reach the equilibrium number A'eq = {Rc / RaY^"^ ■ 

We estimate first the capture rate Rc of WIMPs, using where needed the numerical values 
appropriate for the Sun: The total number of interactions is proportional to the number Ni 
of scatters of type i in the celestial body, to the elastic cross section (Tix, the relative velocity 
Urei and the density of WIMPS nx'- Hence Rc should be of the form Rc oc nx X^i ^if^XiVrei- 
Additionally, a factor / = /(fesc, ^^rei) should count if the scattering was successful, i.e. if 
the WIMP velocity after the scattering is smaller than the escape velocity Vesc from the Sun. 
The function / should be dimensionless, increase for increasing Vesc and can depend only on 
the square of the velocities. Hence / = /(^^esc/'^^rei)' ^ more careful analysis shows that 
indeed / cx fesc/'^^rei- Combining all factors, the capture rate follows as 



Rc = Anx^X.^i^i^ (8.14) 

Vrcl 



V ^ — > 



with A K, 2 for the Sun. To obtain a numerical estimate, we neglect helium and metals 
and set = Mq / mn ■ Using the rotation velocity of the Sun around the Galactic center as 
relative velocity v^ei ~ 220km/h, Vesc = {GMq/ RqY^"^ « 617km/s and nx = p\oc/^x with 
Pioc ~ 0.3GeV/cm^, it follows 

Next we estimate the annihilation rate Ta- Generically, the assumption that the WIMP is 
a thermal relic fixes also the annihilation cross section at present to (cJann^') ~ 3 x lO^^^cm^/s. 
(This could be avoided if s-wave annihilation is suppressed by symmetry reasons.) Since the 
captured WIMPs orbit many times in-between interactions, they can be characterized by a 
global temperature T and their density follows a barometric formula, 

nxi^^) = exp(— mx</'(r)/T) (8.16) 
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determined by the local gravitational potential (j){r). Hence 

Ta = {cT.nnv) 47r / drr^n^xir) . (8.17) 



Determining first Ra = Ta/N'^, we see that this rate is fixed by (cxann^^) and the ratio V2/V1 
of "effective volumes" defined as 

Va=47r/ drr'^exp{-amx(l){r)/T) . (8.18) 
Jo 

To obtain a numerical estimate, we use that WIMPs for mx ^ are concentrated in the 
center of the Sun. Thus we approximate the density by p{r) ~ p(0) ~ 150g/cm^, set T equal 
to the central temperature of the Sun, T = Tc ^ 1.4 x 10'' K, and use Rq — > 00. Then the 
integrals ()8.18p can be performed and one obtains 

y.=6.5xlO-cm'fl55^)''' (8.19) 



a mx 



and thus 

,V2 00 /100GeV\3/2 

RA = {a;nnv)j^ = 6.5xlO^^ . 8.20 

V{ V rnx J 

Now we can determine the equilibration time as r = {RcRaY^/'^ ~ 1.4 X 10^6 s. bmce r is 
much smaller than the age of the Sun, we can use Ta = Rc- Denoting with fy the number 
of neutrinos produced per annihilation, we obtain as our final result for the neutrino flux (^^^ 
from WIMP annihilations in the Sun 

<^ = ^ 5 X 10-8cm-2 s-i sr-i (8.21) 

at the Earth distance d = 1 AU. Comparing this flux with expected neutrino intensities and 
the sensitivity of neutrino telescopes discussed in Sec. 5.5 makes it clear that the detection 
of WIMP annihilations is as challenging as the one of astrophysical neutrinos. The chances 
for success depend strongly on the WIMP mass and its annihilation spectrum. Only if the 
energy of the neutrinos is sufficiently high, Ey > 100 GeV, the direction of the produced muon 
provides a reasonable signature. 



8.4 Summary of possible and suggested signals for DM 

The different annihilation channels and the DM distribution offers several different possible 
signals to detect DM annihilations: 

• All particles from DM annihilations in the Sun or the Earth are absorbed except neu- 
trinos. Apart from a small flux of neutrinos produced by cosmic ray interactions in the 
solar atmosphere, no other neutrinos with energy E > GeV are expected from the Sun. 
The direction provides however only for Ey > 100 GeV a reasonable signature. 

• Cold DM has small velocities and thus XX 27 produces a sharp photon line at 
E^ = mx- This is a smoking-gun for DM, but since the DM particle is most likely 
electrically neutral, the lowest order contributions to these processes are box diagrams, 
and the branching ratio is therefore very small. 
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Figure 8.4: EGRET data on diffuse photon background together photons expected from CR 
interactions (yellow) and a putative DM signal (red) according Ref. [31j . 

• The main production channel for photons are decays of neutral pions. The main problem 
here is to separate a potential DM contribution from diffuse astrophysical Galactic and 
extragalactic backgrounds, cf. Fig. 18.41 More promising might be to look at nearby 
dwarf galaxies that have with a larger DM fraction. 

• Antiprotons could become visible below ~ 1 GeV, but without any specific signatures. 
Thus the identification of DM as source of these antiprotons requires again a precise 
modeling of antiproton flux produced by cosmic rays. 

Various astrophysical observations have been suggested as signals for DM annihilations: 

• The Integral satellite sees a strong positron annihilation line from the Galactic bulge. 
An explanation by DM annihilations into electron-positron pairs, XX e~^e~ , requires 
that the electron-positron pair gets only little kinetic energy. Hence the DM particle 
should be very light, nie < mx ^ lOMeV - which is not excluded but difficult to achieve. 

• The "EGRET excess", i.e. the red region in Fig. 18. 4^ is a possible surplus of diffuse 
gamma-rays compared to the predictions in the simplest models for the propagation of 
cosmic rays. This potential excess has been explained both by DM annihilations in the 
Milky Way and as the sum of DM annihilation in other galaxies. 

• The "WMAP haze" , a potential excess of synchrotron radiation from the Galactic center 
in the WMAP data has been explained as synchrotron radiation of electron produced 
in DM annihilations. 

• HESS, an atmospheric Cherenkov telescope, has observed TeV 7-rays from the Galactic 
center. The flux extends however power-law like up to the sensitivity limit (~ 30 TeV) 
and is thus difficult to combine with the upper limit on M for a thermal relic. 

For references to the original literature and a brief critical review see Ref. [S2\ . 
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Congratulation to the readers who worked through these notes until the end. Certainly many 
questions have been left unanswered and some new ones have popped up. Still I hope that 
these notes could provide a first introduction into high energy cosmic ray physics. For further 
studies I recommend the following two excellent references, 

[El] V. S. Berezinsky, S. V. Bulanov, V. A. Dogiel, V. L. Ginzburg, V. S. Ptuskin, Astro- 
physics of Cosmic Rays (North-Holland, Amsterdam 1990), 

[E2] T. K. Gaisser, Cosmic Rays and Particle Physics (Cambridge University Press, Cam- 
bridge 1990). 

As the titles indicate, the focus of the first book is on astrophysical aspects of cosmic rays, 
while the second discusses additionally the interactions of cosmic rays in the atmosphere and 
the formation of extensive air showers. The presentation of several topics I have given was 
heavily influenced by these books. I omitted generally references to the original literature, 
apart from references to the sources of the reproduced figures, and refer also in this respect 
to the books of Berezinsky et al. and Gaisser. Additionally, the following review articles may 
be helpful: 

1. V. Berezinsky, "Ultra high energy neutrino astronomy," Nucl. Phys. Proc. Suppl. 151, 
260 (2006) |ar Xiv:astro- ph/0505220|. 

2. P. Lipari, "Perspectives of high energy neutrino astronomy," Nucl. lustrum. Meth. A 
567, 405 (2006) [arXiv:astro-ph/0605535| . 

3. V. Berezinsky, "Transition from galactic to extragalactic cosmic rays," to appear in the 
proceedings of 30th ICRC 2007, .arXiv:0710.2750., [astro-ph]. 

4. M. Kachelriefi, "Status of particle physics solutions to the UHECR puzzle," Comptes 
Rendus Physique 5, 441 (2004) |arXiv:hep-ph/0406174 |. 

5. G. Jungman, M. Kamionkowski and K. Griest, "Supersymmetric dark matter," Phys. 
Rept. 267, 195 (1996) |arXiv:hep-ph /9506380|. 

It is a pleasure to thank all my collaborators in the field of cosmic ray physics, but in 
particular Roberto Aloisio, Venya Berezinsky, Sergey Ostapchenko, Dima Semikoz, Pasquale 
Serpico, and Ricard Tomas, as well as Giinter Sigl for many illuminating discussions. I am 
grateful to Sergey Ostapchenko and Pasquale Serpico for reading (parts of) the manuscript 
and pointing out several errors and misunderstandings contained in the text. 

Last but not least, I would like to thank Kimmo Kainulainen for inviting me to lecture at 
the Jyvasklya Summer School and for the excellent working conditions provided there. 



79 



Bibliography 



[1] E. S. Seo et al, "Measurement of cosmic-ray proton and helium spectra during the 1987 
solar minimum," Astrophys. J. 378, 763 (1991); M. Nagano et al, "Energy spectrum of 
primary cosmic rays between 10^^'^ and 10^^ eV," J. Phys. GIO, 1295 (1984); M. Nagano 
et al., "Energy spectrum of primary cosmic rays above 10^^ eV determined from extensive 
air shower experiments at Akeno," J. Phys. G18, 423 (1992); N. L. Grigorov et al, 
"Energy spectrum of primary cosmic rays in the 10^^ — 10^^ eV according to the data of 
Proton-4 measurements," Proceedings 12th ICRC, 1, 1760 (1971). 

[2] V. Hess, "Uber die Beobachtung der durchdringenden Strahlung bei sieben Freiballon- 
fahrten," Phys. Z. 13, 1084 (1912). 

[3] P. Desgrolard, M. Giffon, E. Martynov and E. Predazzi, "Exchange-degenerate Regge 
trajectories: A fresh look from resonance and forward scattering regions," Eur. Phys. J. 
C 18, 555 (2001) | arXiv:hep-ph/0006244] . 

[4] W.-M. Yao et al, J. Phys. G 33, 1 (2006). 

[5] M. H. Israel et al. [ACE collaboration], "Isotopic Composition of Cosmic Rays: Results 
from the Cosmic Ray Isotope Spectrometer on the ACE Spacecraft," Nucl. Phys. A758, 
201 (2005). 

[6] M. E. Wiedenbeck et al. [ACE collaboration], "The level of solar modulation of Galactic 
cosmic rays from 1997 to 2005 as derived from ACE measurements of elemental energy 
spectra," Proc. "19th International Cosmic Ray Conference" , Pune, India, 2005, 2, p. 
277 (2005). 

[7] T. Antoni et al. [The KASCADE Collaboration], "Large scale cosmic-ray anisotropy with 
KASCADE," Astrophys. J. 604, 687 (2004) [arXiv: astro-ph'/Cm 2375| . 

[8] Figure from |http : / /hep . bu . edu/~superk/ ew- ef f ect . htinl[ 

[9] K. H. Kampert, "Cosmic rays from the knee to the ankle: Status and prospects," Nucl. 
Phys. Proc. Suppl. 165, 294 (2007) [arXiv:astro-ph/0611884]. 

[10] V. S. Berezinsky, S. I. Grigorieva and B. I. Hnatyk, "Extragalactic UHE proton spectrum 
and prediction for iron-nuclei flux at 10**8-GeV to 10**9-GeV," Astropart. Phys. 21, 
617 (2004) iajXiv:astro-ph /0403477] . 

[11] F. A. Aharonian, W. Hofmann, A. K. Konopelko, and H. J. Volk, "The potential of the 
ground based arrays of imaging atmospheric Cherenkov telescopes," Astropart. Phys. 6, 
369 (1997). 

[12] R. D. Blandford, "Acceleration of ultrahigh-energy cosmic rays," Phys. Scripta T85, 191 
(2000) [arXiv:astro-ph/ 9906026|. 



80 



Bibliography 



[13] P. Blasi, "On the origin of very high energy cosmic rays," Mod. Phys. Lett. A 20, 3055 
(2005) |arXiv:astro-ph/0511235|; M. Ostrowski, "Cosmic ray acceleration at relativistic 
shocks, shear layers,. . . ," to appear in New Astron. Rev., arXiv:0801.1 339 [astro-ph]. 

[14] F. A. Aharonian et al. [The HESS Collaboration], "High-energy particle acceleration in 
the shell of a supernova remnant," Nature 432, 75 (2004) arXiv:astro-ph/0411533j . 

[15] Figure from William Keel, University of Alabama, 

|http : //www . astr . ua . edu/keel/agn/mkn421 . html 

[16] B. Rossi and K. Greisen, "Cosmic-ray theory," Rev. Mod. Phys. 13, 240 (1941). 

[17] A. Miicke et al, "Monte Carlo simulations of photohadronic processes in astrophysics," 
Comput. Phys. Commun. 124, 290 (2000) [astro-ph /99b3478]. 

[18] M. Kachelriefi, P. D. Serpico and M. Teshima, "The galactic magnetic field as 
spectrograph for ultra-high energy cosmic rays," Astropart. Phys. 26, 378 (2006) 
[arXiv:astro-ph/0510444j . 

[19] M. KachelrieB and D. V. Semikoz, "Reconciling the ultra-high energy cosmic 
ray spectrum with Fermi shock acceleration," Phys. Lett. B 634, 143 (2006) 
[arXiv:astro-ph/05101 88| . 

[20] V. Berezinsky, A. Z. Gazizov and S. L Grigorieva, "On astrophysical solution to ultra 
high energy cosmic rays," Phys. Rev. D 74, 043005 (2006) [arXiv:hep-ph/0204357 |. 

[21] A. Cuoco et al., "The footprint of large scale cosmic structure on the ultrahigh energy 
cosmic ray distribution," JCAP 0601, 009 (2006) [ arXiv:astro-ph/0510765 1. 

[22] M. Kachelriefi and D. V. Semikoz, "Clustering of ultra-high energy cosmic ray arrival 
directions on medium scales," Astropart. Phys. 26, 10 (2006) [arXiv: astro-ph /05 1 2498 J . 

[23] S. Mollerach et al. [PAO Collaboration], "Studies of clustering in the arrival directions of 
cosmic rays detected at the Pierre Auger Observatory above 10 EeV," to appear in Proc. 
"30th International Cosmic Ray Conference" , Merida, Mexico, 2007, larXiv:0706. 17491 
[astro-ph] . 

[24] A. Cuoco et al, 'Clustering properties of ultrahigh energy cosmic rays and the search for 
their astrophysical sources," to appear in the Astrophys, J.. larXiv:0709.2712l [astro-ph]. 

[25] J. Abraham et al. [PAO Collaboration], "Correlation of the highest energy cosmic rays 
with nearby extragalactic objects," Science 318, 939 (2007). larXiv:0711. 22561 [astro-ph]. 

[26] M. Kachelriefi and R. Tomas, "High energy neutrino yields from astrophysical sources. L 
Weakly magnetized sources," Phys. Rev. D 74, 063009 (2006) [arXiv:astro-ph/0 606406 1 . 

[27] K. Mannheim, R. J. Protheroe and J. P. Rachen, "On the cosmic ray bound 
for models of extragalactic neutrino production," Phys. Rev. D 63, 023003 (2001) 
[arXiv:astro-ph/9812398i . 

[28] R. Gandhi, C. Quigg, M. H. Reno and L Sarcevic, "Neutrino interactions at ultrahigh 
energies," Phys. Rev. D 58, 093009 (1998) ^arXiv:hep-ph/9807 264 1 . 



81 



Bibliography 



[29] H. J. Drescher, A. Dumitru and M. Strikman, "High-density QCD and cosmic ray air 
showers," Phys. Rev. Lett. 94, 231801 (2005) |arXiv:hep-p h/0408073]. 

[30] K. Griest and D. Seckel, "Cosmic Asymmetry, Neutrinos and the Sun," Nucl. Phys. B 
283, 681 (1987) [Erratum-ibid. B 296, 1034 (1988)]. 

[31] W. de Boer, "Indirect evidence for WIMP annihilation from diffuse galactic gamma 
rays," arX iv:astro-ph/0412620 

[32] D. Hooper, "Indirect Searches For Dark Matter: Signals, Hints and Otherwise," to appear 
in the proceedings of SUSY07. !aFXiv:0710. 2062 [hep-ph]. 



82 



